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1. Introduction 

Quantum chroniodynamics, the theory of the strong interactions, is a non- 
AbeUan gauge theory based on the gauge group SU (3). It was first pointed 
out by 't Hooft [1,2] that many features of QCD can be understood by 
studying a gauge theory based on the gauge group SU{N) in the hmit 
N ^ oo. One might think that letting N ^ oo would make the anal- 
ysis more complicated because of the larger gauge group and consequent 
increase in the number of dynamical degrees of freedom. One might also 
think that SU{N) gauge theory has very little to do with QCD because 
A'' = DO is not close to A = 3. However, we will soon see that SU (N) gauge 
theory simplifies in the N ^ oo limit, that the true expansion parameter 
is 1/N, not N, and that the 1/N expansion is equivalent to a semiclassi- 
cal expansion for an effective theory of color singlet mesons and baryons. 
Results for QCD can be obtained from the N oc' limit by expanding in 
1/N = 1/3, and are in good agreement with experiment. 

To decide whether 1/N is a small expansion parameter for QCD requires 
further analysis. In QED, as Wittcn has remarked, the coupling constant 
e = V^Tra = 0.30, which is not very different from 1/3. Anyone who 
has actually computed radiative corrections in QED knows that the true 
expansion parameter is not e, but is closer to a/An w 10"'^, which is 
much smaller than e. By the end of these lectures you will see several 
examples which show that the expansion parameter for QCD is 1/N = 1/3. 
While not as small as the QED expansion parameter a/Air, 1/N is still a 
useful expansion parameter for QCD. 1/A corrections are comparable in 
size to flavor SU{3) breaking corrections due to the strange quark mass, 
and expanding in flavor SU (3) breaking is well-known to be an extremely 
useful expansion in QCD. Furthermore, we will find many examples where 
the 1/N term vanishes, so that the first correction is of order 1/N^. In 
such cases, one can make predictions at the 10% level. This is a level of 
computational accuracy in low-energy hadronic physics that is difficult to 
match using other techniques. 

In these lectures, I will concentrate on the large N expansion for QCD, 
and in particular, on trying to obtain QCD results that can be compared 
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with experimental data. Sections 2 3 and sections 4.6, 5 5.1 of these lec- 
tures are based on the treatments by Coleman [3], and Witten [4,5], re- 
spectively. Large N expansions have also been used to study other field 
theories, such as the 0(N) model, CP'^ model, etc. They provide insight 
into quantum field theory dynamics, and have many applications in high 
energy physics and statistical mechanics. They have been extensively used 
in recent years to study matrix models. These topics have been discussed 
in previous Les Houches summer schools, and will not be repeated here. 
A good reference on large A'' methods is the compilation by Brezin and 
Wadia [6]. 

2. The Gross-Neveu Model 

The Gross-Neveu model [7] is an interesting 1 + 1 dimensional field theory 
that can be studied using the 1/A'' expansion. The model is asymptotically 
free with a spontaneously broken chiral symmetry, and so shares some 
dynamical features with QCD. It will provide a useful warm-up exercise 
before we tackle the much more difficult problem of large N QCD. 
The Gross-Neveu Lagrangian is 

L = V^i^V+^(?^V')', (2.1) 

where 7/;°, a = 1, . . . , A arc TV Dirac fields, and a sum on N is implicit in 
the notation, so that tpip = ^°-ip°-, etc. In 1 -|- 1 dimensions, Dirac fields 
are two-component spinors, and have mass dimension 1/2. A is a dimen- 
sionless coupling constant. Equation (2.1) is invariant under an SU{N) 
flavor symmetry on the ip^s, 

where U is an SU{N) matrix, and also invariant under a discrete chiral 
symmetry 

^ — > ip — > —'075; V'V' — * —tpip. (2-2) 

Equation (2.1) is the most general possible Lagrangian invariant under 
these symmetries with terms of dimension less than or equal to two, and 
so describes a renormalizablc field theory in 1 + 1 dimensions. The discrete 
chiral symmetry eq. (2.2) forbids a mass term, so the fermions are massless 
at any finite order in perturbation theory, and no mass counterterm is 
needed to regulate the ultraviolet divergences. 

The basic interaction vertex is shown in fig. 1. Consider the process 
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a a 




b b 



Fig. 1. The four-Fermi vertex of the Gross-Neveu model, o, b are flavor labels. 



a b a b 




a b a b a b 

(c) (d) (e) 



Fig. 2. Low-order diagrams for the scattering amplitude a + a ^ b + b, with a ^ b. 
Virtual flavors c and d are summed over. 



a -|- o — > h + with a ^ h. Some graphs contributing to this scattering 
ampHtude arc shown in fig. 2. The leading order diagram fig. 2(a) is of 
order A. The one-loop correction fig. 2(b) is of order A^. The intermediate 
fermion flavor c in the one-loop correction fig. 2(c) is arbitrary and must 
be summed, so the graph is order A^A'^. Similarly, figs. 2(d,c) arc of order 
X^N"^. One can clearly see that the perturbation series does not have a 
well-defined limit as A?' — > oo, because the radiative corrections grow with 
powers of N . 

One can obtain a well-defined large N limit by rescaling the coupling 



8 



A. V. Manohar 



constant A. Define A = g'^ /N, and take the limit N ^ oo with g fixed, 
so that A is of order 1/N. The graphs considered in fig. 2 then give a 
sensible expansion for the scattering amplitude — the lowest order term 
is g'^/N, and the correction terms arc g'^/N'^, g'^/N, g^/N and g^/N for 
figs. 2(b-e), respectively. This shows that the perturbation series gives 
a scattering amplitude of the form {1/N)f{g^,l/N), and has the large N 
limit (l/iV)/(.g2, 0). In the large N limit, dia grams figs. 2(a,c c) contribute 
to f{g^,0), but diagram fig. 2(b) is omitted. There is some simplification 
of the diagrammatic expansion as N ^ oo, but the limit is still highly 
non-trivial. 

The Gross-Neveu Lagrangian is 

2 

L = ^i$^+-^(^i,)\ (2.3) 

when written in terms of g. One way to understand the power of 1 /N in the 

interaction term is to note that iljijj / \fN produces a flavor singlet il)ip state 
with unit amplitude, since there is an implicit sum over N flavors in il)tp. 
This is like in quantum mechanics, where a state which is the sum of N 

orthonormal states with equal amplitude, = a + \2) + . . . + \N)) , 

has normalization constant N \ a\ = 1 to have unit norm. Then (iptp) 
should have a coeSicient of order 1/N, so that tl^tp scattering in the flavor 
singlet channel has an amplitude of order unity. 

One can now study the perturbation series in g for the Lagrangian 
eq. (2.3) in the N ^ oo limit. The diagrammatic expansion for the Gross- 
Neveu model simplifies in the large N limit. For example, wc have seen 
that fig. 2(b) can be neglected. The simplifications are sufficient to allow 
one to obtain exact results, though this might not yet be apparent. To 
make the large N analysis more transparent, it is convenient to introduce 
an auxiliary field a, and write the Lagrangian eq. (2.3) as 

L=lbiStb + (7lbib- ^a"^. (2.4) 
2g^ 

The Lagrangian is quadratic in a, so integrating over a is equivalent to 
minimizing the Lagrangian with respect to a, which gives 

a = j^ip-ip- 

Substituting the answer back in eq. (2.4) gives the original form of the 
Lagrangian, eq. (2.3). 

The analysis of the large N limit is simpler using the modified La- 
grangian eq. (2.4). The Feynman rules are given in fig. 3, and the scattering 
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Fig. 3. Feynman rules for the Gross-Neveu Lagrangian eq. (2.4). The V field is repre- 
sented by a solid line, and the auxiliary field cr by a dashed line. 




(a) (b) (c) 

Fig. 4. The graphs of fig. 2(a,c) using the auxiliary field form of the Lagrangian. 



graphs of fig. 2 are now represented as in fig. 4. Tlie auxiliary field rep- 
resentation allows one to obtain the iV-counting rules for diagrams. Con- 
sider the graphs of fig. 4 with all external fermion lines removed. The 
resulting graphs are generated by an effective action Les{<^) which con- 
tains only external a lines. This effective action is obtained by evaluating 
the fermion functional integral using the Lagrangian eq. (2.4). The La- 
grangian is quadratic in the fermion fields, so Leg{a) is given exactly by 
the sum of diagrams in fig. 5. The first term is the tree- level inverse prop- 
agator —N/2g'^a^, and the remaining terms are the one-loop corrections. 
Each diagram in fig. 5 is of order N — the one-loop terms have N fermions 
in a closed loop, and the tree-level term is explicitly of order A''. Thus the 
effective action can be written as 

L^s{<7,g,N)=NL^e{(T,g). (2.5) 

It is now straightforward to determine the power of A'^ in any Feynman 
graph with only external a lines. Each term in the Lagrangian eq. (2.5) is 
of order N. Each interaction vertex has a factor of N, and each propagator 
has a factor of 1/A'^, since a vertex is a term in the Lagrangian, and a 
propagator is the inverse of the quadratic terms in the Lagrangian. Thus a 
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diagram is proportional to jy'^-i-E ^ where V is the number of vertices, / 
is the number of internal a lines, and E is the number of external a lines. 
Factors of are included for external a propagators since the physical 
scattering amplitudes have only external fermion fields, and all the a lines 
are actually internal a lines in the full diagram including tp propagators. 
An example is shown in fig. 6. Diagrams (a) and (b) can be redrawn as (c) 
and (d), where the blob is an interaction vertex in NLcs- The A'^-counting 
formula then shows that fig. 6(c) is of order N^^'^^'^ = 1/^, and fig. 6(d) 
is of order A^^"^"^ = These are also the A^-counting rules for the 

original diagrams figs. 6(a) and (b), respectively. For any Feynman graph, 
one has the identity 

V-I + L = l, (2.6) 

so that a Feynman diagram is of order 

N^-L-E_ (2.7) 

The minus signs in front of L and E in eq. (2.7) arc important, since they 
imply that additional loops or external lines bring an additional suppression 
of powers of 1/iV, rather than enhancements by powers of N. This proves 
that the theory has a sensible 1/N expansion. Since the minimum number 
of external a lines is at least two, the maximum power of A'' is —1. 

The effective Lagrangian Leff can be computed from the bubble sum in 
fig. 5. It will be computed here only in the limit where all a lines carry zero 
momentum, where it reduces to the effective potential V{cr). The effective 
potential is given by the sum of diagrams in fig. 5, 
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(a) 



(b) 



(c) 



(d) 



Fig. 6. Examples of AT-counting. Graphs (a) and (b) are equivalent to (c) and (d), where 
the blobs are interaction vertices in NLgg. 



There are only even powers of a because of the discrete symmetry a —a. 
The first term in eq. (2.8) is the tree-level amplitude. The second term is the 
sum of loop graphs. The —A'' is from the flavors of fermions in the loop, 
l/2r is the symmetry factor for the graph, and the term in parentheses is 
the product of the fermion propagator ip/p^ and the vertex ia. Performing 
the trace, using the identity 



2r 



r—l 



1 



log (1 - ar^) 



and analytically continuing to Euclidean space gives 



V 



N 
2? 



N 



log 



1 + ^ 



(2.9) 



Regulating the loop integral using dimensional regularization in the MS 
scheme gives 



V = N 



2ff2 



a 
4tt 



1 



(2.10) 



The eff'ective potential V{a) satisfies the renormalization group equation 

V{a)=0 (2.11) 



'9 ^, , d 
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Substituting eq. (2.10) into eq. (2.11) gives 

7.(g) = 0, /3(.9) = -|^- (2.12) 

These are the exact anomalous dimension and /3-function to all orders in 
g in the N oo Umit. The Gross-Neveu model is an asymptotically free 
theory, since the /3-function is negative. 

The Gross-Neveu model also exhibits spontaneous symmetry breaking. 
The extrema of the effective potential eq. (2.10) are at 

(7 = 0, a = zb^e-'^/s' = ±ao, 
at which V has the values 

2 

V{0)=0, V{±ao) = -N^ <0, 

47r 

so that the global minima of the potential are a = ±ao- The discrete 
symmetry eq. (2.2) is spontaneously broken, since 

and the two minima icjo are mapped into each other under this broken 
symmetry. The fermions get a mass m = ao, since the Yukawa coupling is 
atjjip, and there is no wavefunction renormalization of the a field. 

The key simplification of the large N limit was that the diagrams of 
the theory reduced to a subset, fig. 5, which could be summed exactly to 
give an effective Lagrangian L^s- The large A'' limit is the same as the 
semiclassical limit for Lgff{a). This is evident from the overall factor of N 
in Lef[{a), eq. (2.5). The form of the functional integral 

J Va e'^^«/'' = / e^^^««/'^ 

shows that an expansion in h is equivalent to an expansion in 1/A''. 
The Gross-Neveu Lagrangian eq. (2.3) can be written as 

L = N'^i^^ + Ng'^ C^^f , (2.13) 

using rescalcd fcrmion fields tp = y^^*- This also has an overall factor 
oi N, so one might naively think that the large N limit is the same as 
the semiclassical limit for the Lagrangian (2.13). This is incorrect, because 

the terms in the Lagrangian have hidden N dependence, because there are 
N flavors of and Feynman diagrams have factors of N from the flavor 



Large N QCD 



13 



index sums. The effective action S'ofr (cr) has no hidden factors of N, since 
(7 is a single component flavor singlet field. In this case, the overall factor 
of A'' does imply that the large A'' and semiclassical limits are the same. 
The effective action Scs (cr) for the composite field a is obtained by adding 
tree and loop graphs in the original Gross-Neveu theory, and so contains 
quantum corrections in the Gross-Neveu Lagrangian eq. (2.3). The large 
N limit of the Gross-Neveu model is thus equivalent to the semiclassical 
expansion of an effective theory of flavor singlet a fields ( "mesons" ) . The 
a effective action includes quantum corrections in the Gross-Neveu model, 
so the large A'^ limit is not the same as the semiclassical limit of the original 
Gross-Neveu model. A similar result holds for QCD. We will see that the 
large A'' limit of QCD is the same as the semiclassical limit of an effective 
theory of color singlet mesons and baryons. 

Problem 2.1 (Unitarity Bound) 

Show that the a + a ^ b + b amphtude must be of order 1/N (or smaller) 
to avoid violating unitarity in the large N limit. 

Problem 2.2 

Prove eq. (2.6). 

Problem 2.3 (EfTective Potential in the MS Scheme) 

Evaluate the effective potential eq. (2.9) by analytically continuing the 
momentum integral to 2 — 2e dimensions. You can apply the familiar rules 
for Feynman graphs in D dimensions by using the identity 

da V PV °^ V ' 
Problem 2.4 (1/N Corrections to V) 

The effective potential for a contains two loop corrections, such as fig. 7(a). 
In the method of calculation outlined above, where one first computes the 
fermion functional integral, these terms are obtained by computing loop 
graphs using Loir- a^s in fig. 7(h). These are suppressed by 1/N, and were 
neglected above. Include these (unknown) order g'^/N terms in V, and 
repeat the derivation of eq. (2.12). Assume that 7<j (g) starts at order g^, 
and (3 (g) starts at order g^. Show that one obtains 

7a{9) = ^0{g'), 
H9) = -^-97A9) + io{g'). 
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Fig. 7. A two loop correction to the cfTcctivc potential V{a). The correction in the theory 
with both i/) and cr fields is given by (a). The same diagram in the theory obtained by 
integrating out tp is given by (b). 



3. QCD 

3.1. N- Counting Rules for Diagrams 

The analysis of the AT-counting rules for QCD is more complicated than 
that for the Gross-Neveu model studied in the previous section. The main 
reason for this is that gluons transform under the adjoint representation of 
the gauge group, rather than the fundamental representation. In the Gross- 
Neveu model, the dynamics could be rewritten in terms of a singlet field 
a = tjj'4'- In QCD, one can construct an infinite number of gauge singlets, 
e.g. TrF^^, TrF^^^, . . ., TrF^, from the gluon field-strength tensor F^^. 

The theory we will study is an SU{N) gauge theory with Np flavors 
of fermions (quarks) in the fundamental representation of SU{N). The 
gauge field is an N x N traceless hermitian matrix, = A^T"^, and the 
covariant derivative is 

The matrices T"* are normalized so that 
2 

The coupling constant has been chosen to be g/^/N, rather than g, because 
this will lead to a theory with a sensible (and non- trivial) large N limit. 
The field strength is 



F^v — dfxAi, di/A^ -\- i—j= \A^, Ad] , 



Large N QCD 



15 



and the Lagrangian is 

Nf 

TV F.^.F^" A- 

2 



1 Nf ^ 

L = --TVF^,i^''^ + ^V^fe(ip-mfc)Vfe. (3.1) 

The large iV Umit will be taken with the number of flavors Np fixed. It 
is also possible to consider other limits, such as N ^ oo with Np/N held 
fixed [8]. 

One way to understand the g/VN scaling of the coupling constant is to 
look at the QCD /3-function, 

using the conventionally normalized coupling constant. This equation does 
not have a sensible large A'' limit since bo is order A''. Replacing g by g/VN 
in eq. (3.2) gives 

dg (II 2Np\ g^ , 

The /3-function equation now has a well-defined limit as iV od. The 
Np term is suppressed by and we will soon see that all fermion loop 
effects are suppressed. The scale parameter of the strong interactions, 
AqcD) is held fixed as ^ oo, since N drops out of the equation for the 
running of g. Thus the large N limit for QCD with the coupling constant 
scaling like 1/Vn is equivalent to taking the limit N ^ oo holding the 
string tension, or a meson mass such as the p mass, fixed. 

To analyze the A^-counting rules for QCD, one needs a simple way to 
count the powers of A^ in a given Feynman diagram. This can be done with 
the help of a trick due originally to 't Hooft. The quark propagator is 

r{x)W{y)) = S'''S{x-y). (3.3) 

This is represented diagrammatically by a single line (fig. 8(a)), and the 
color at the beginning of the line is the same as at the end of the line, 
because of the (5"^ in eq. (3.3). The gluon propagator is 

{A^ {x)A^ {y))=5^''D^, {x - y) , 



whore A and B arc indices in the adjoint representation. Instead of treating 
a gluon as a field with a single adjoint index, it is preferable to treat it 
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(a) (b) 




(c) (d) (e) 



Fig. 8. 't Hooft's double line notation. The lower diagram shows each QCD propagator 
or interaction vertex in the double line notation. 



as an iV X A'' matrix with two indices in the N and N representations, 
b = A^ {^^) b- The gluon propagator can be rewritten as 

(A-b (x) AU {y)) = D^^ {x - y) - , 

where the identity 

{T^)\ {T^)\ = - 5'i {SU{N)) (3.4) 

has been used. The corresponding identity for U{N) is 

(T^)\M% = i<5S<5,^ (C/(7V)), (3.5) 

where the U{1) generator has the same normalization as the SU{N) gen- 
erators. It is convenient to use the U{N) identity eq. (3.5) instead of the 
SU{N) identity eq. (3.4) for analyzing the A'^-dependence of Feynman dia- 
grams. The correct SU (N) propagator is given by including an additional 
U{1) ghost field that cancels the extra U{1) gauge boson in U{N). The 
effects of the U{1) gauge boson are l/N"^ suppressed, as we will see later. 
In most applications, the difference between U{N) and SU{N) will turn 
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Fig. 9. Color flow for a four-gluon interaction which has two color traces. 



out to be unimportant. The reason is that, typically, what one can prove 

is that a certain amplitude first occurs at some order in 1/A''. For such 
computations, the numerical size of the 1/N corrections is irrelevant. 

The U{N) gluon propagator can then be represented using 't Hooft's 
double line notation, as in fig. 8(b). The color lines represent the N and N 
indices a and b on A^^) and color is conserved during propagation because 
of the i5-function structure in eq. (3.5). The gauge- fermion vertex ip 
is shown in fig. 8(c). The double line notation provides a simple way to 
keep track of the color index contractions in a Feynman graph. 

The three and four gauge boson vertices arise from the Tr F^vF^^'^ gluon 
kinetic energy. Each kinetic energy term is a single trace over color. The 
three-gluon vertex arises from terms such as 

where the color indices have been shown explicitly. It is represented using 
the double line notation as fig. 8(d), and the four gluon vertex as in fig. 8(e). 
It is important that all the interactions arise from a single color trace — 
otherwise one could have color flow in a four-gluon vertex as in fig. 9, where 
the diagram can be broken up into two disconnected color flows. 

Every Feynman graph in the original theory can then be written as a 
sum of double line graphs. Each double line graph gives a particular color 
index contraction of the original diagram. An example of a Feynman graph 
and one of its double line partners is shown in fig. 10. One Feynman graph 
can give rise to several double line graphs. For example, the three-gluon 
vertex is given by eq. (3.6), plus a term with <-> A^. The complete 
three-gluon vertex in the double line notation is represented in fig. 11. 

One can think of each double line graph as a surface obtained by gluing 
polygons together at the double lines. Since each line has an arrow on it, 
and double lines have oppositely directed arrow, one can only construct 
orientable polygons in an SU{N) theory. For SO{N), the fundamental 
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Fig. 10. A Fcynman graph, and one of its partner double line graphs, representing a 
particular color index contraction. This is an example of a planar diagram. 




representation is a real representation, and the lines do not have arrows. 
In this case, it is possible to construct non-orientable surfaces such as Klein 
bottles. 

To compute the A/'-dependence requires counting powers of N from sums 
over closed color index loops, as well as factors of 1/ \/N from the explicit N 
dependence in the coupling constants. It is convenient to use a rescaled La- 
grangian to simplify the derivation of the A/'-counting rules. Define rescaled 
gauge fields gA/ ^/N ^ Aso that the covariant derivative is D/^ = df^^+iAf^, 
and rescaled fermion fields tp \fN tp so that the Lagrangian becomes 

L = N 

The Lagrangian has an overall N, but the theory does not reduce to a 
classical theory of quarks and gluons in the N ^ oo limit, because the 
number of components of ip and grows with N. 

One can read off the powers of TV in any Feynman graph from eq. (3.7). 
Every vertex has a factor of N, and every propagator has a factor of 1/A^. 
In addition, every color index loop gives a factor of N, since it represents 
a sum over N colors. In the double line notation where Feynman graphs 
correspond to polygons glued to form surfaces, each color index loop is the 



Nf 



fe=i 



(3.7) 
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edge of a polygon, and is the face of the surface. Thus one finds that a 
connected vacuum graph (i.e. with no external legs) is of order 



where V is the number of vertices, E is the number of edges, F is the 
number of faces, and x = V — E + F \s & topological invariant known as 
the Euler character. For a connected orientable surface 



where h is the number of handles, and h is the number of boundaries (or 
holes). For a sphere /i = 0, 6 = 0, x = 2; for a torus, /i = 1, 6 = 0, x = 0. 

A quark is represented by a single line, and so a closed quark loop 
is a boundary. Thus every closed quark loop brings a 1/N suppression. 
The maximum power of N is two, from graphs with h = h = Q. These 
are connected graphs with no closed quark loops, with the topology of a 
sphere. Remove one polygon from the sphere, so that one obtains a sphere 
with one hole. This can be flattened into a diagram drawn on a flat sheet 
of paper, with the hole as the outermost edge. One can then glue back 
the removed polygon by thinking of it as the paper exterior to the diagram 
with infinity identified. Thus the order N"^ graphs are planar diagrams 
with only gluons, that is, they can be drawn on the surface of a sheet of 
paper without having a gluon "jump" over another. All points where gluon 
lines cross have to be interaction vertices. 

A diagram with c connected pieces can be of order N'^'^, and so grows 
with the number of connected pieces. This is not surprising; the sum of 
all diagrams is the exponential of the connected diagrams. The connected 
diagrams are of order N'^, corresponding to a vacuum energy of order N'^, 
which is to be expected since there are O (iV^) gluon degrees of freedom. 
Expanding the exponential gives arbitrary high powers of N . From now 
on, we will restrict A''-counting to the connected diagrams. One can obtain 
the N dependence of a disconnected diagram by multiplying together the 
N dependence of all the connected pieces. 

We will often be interested in correlation functions that depend on prop- 
erties of the quarks, such as masses. The leading graphs that depend on 
quarks must have at least one quark line, and are order N , with h = and 
6=1. One might expect the quark contribution to the vacuum energy to 
be of order TV, since there are N quarks of each flavor. The order A'^ quark 
diagrams have the topology of a sphere with one hole, with the quark loop 
forming the edge of the hole. One can then flatten them out into a planar 
diagram as for gluons. In this case, the order N diagrams are written as 
planar diagrams with a single quark loop which forms the outermost edge 



j^V-E+F ^ j^x^ 



(3.8) 



X = 2 - 2/1 - 6, 



(3.9) 
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Fig. 12. An example of a non-planar diagram. Diagram (a) can be redrawn as (b), but 
it still non-planar since the quark loop must form the outermost edge of the diagram. 
The double line version of the graph is shown in (c) . 



of the diagram. An example of a planar quark diagram is fig. 10. Figure 10 
is of order N, since h = and 6=1. This can also be seen using the 
original A^-counting rules of the Lagrangian eq. (3.1). Each vertex has a 
factor g/^/N, and each closed index loop brings a factor of A^, so the graph 
is of order (l/\/]V)* x = N. Figure 12(a) is not a planar diagram, even 
though it can be drawn as fig. 12(b), because the diagram must be planar 
when drawn with the quark line as the outermost edge. Figure 12(a) is of 
order 1/A', since the vertex factors give (l/y/N)'^ and the color index sum 
gives a factor of A^. Note that for a given number of quark loops (bound- 
aries 6), the expansion is in powers of l/N^, rather than l/N, because of 
the —2 in front of h in eq. (3.9). 

Large N diagrams for QCD look like two-dimensional surfaces. For ex- 
ample, the leading diagram in the pure-glue sector has the topology of a 
sphere, and the leading diagram in the quark sector is a surface with the 
quark as the outermost edge. One can imagine all possible planar gluon 
exchanges as filling out the surface into a two-dimensional world-sheet. It 
has been conjectured that this is the way in which large N QCD might be 
connected with string theory, with planar diagrams representing the lead- 
ing order string theory diagrams. The topological counting rules for the 
1/A' suppression factors in QCD are the same as that for the string cou- 
pling constant in the string loop expansion. The connection between large 
N QCD and string theory has never been made precise. Two major ob- 
stacles are that QCD is neither supersymmetric nor conformally invariant. 
One result in this direction is that the partition ftmction for SU{N) Yang- 
Mills theory (i.e. no quarks) in two dimensions was shown to agree with 
the partition function of a string theory [9] by explicit calculation. The 
connection was possible because Yang-Mills theory in two-dimensions is a 
free field theory, and the partition function only depends on the topology 
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Fig. 13. A connected diagram with one U (1) ghost gauge boson exchange. 



of the background spacetime. To reproduce the correct iV-factors, it was 
necessary to use a modified string theory, in which folds are suppressed. 



3.1.1. U{1) Ghosts 

The corrections due to the difference between using SU{N) and U{N) are 
straightforward to analyze. As mentioned earlier, the SU {N) theory can 
be thought of as a U{N) theory with an additional U{1) ghost gauge field 
to cancel the extra ^7(1) gauge field in U{N). The U{1) ghost field does not 
couple to the U{N) gauge bosons, since the U{1) generator commutes with 
all the U{N) generators. We only need to consider exchanges of the U{1) 
ghost field between quark lines. The additional powers of N due to the U (1) 
ghost are most simply counted using the original form of the Lagrangian 
eq. (3.1), rather than the rescaled version eq. (3.7). Consider a connected 
diagram, with some gluon lines and ghost U{1) lines. The C/(l) ghost does 
not change the color structure of the diagram, so the N dependence from 
the color index loops, etc. is obtained using the counting rules discussed 
above, for the diagram with the C/(l) ghosts erased. In addition, one gets a 
factor of 1/N from each U{1) ghost propagator (see eq. (3.4)), and a factor 
ofl/N from the two coupling constants at the two ends of each gauge boson 
line. Thus each C/(l) ghost brings a suppression. The only subtlety 

in this argument is that U{1) exchange can turn an otherwise disconnected 
diagram into a connected diagram. Thus the leading diagram with one 
f/(l) ghost has two quark loops, as in fig. 13, and is order N x N x 1/iV^ = 
O (1). This is only 1/A^ suppressed relative to the leading connected quark 
diagrams, which are order N. However, the effect of the U{1) ghost in 
fig. 13 is to precisely cancel the corresponding graph with a U{N) boson, 
since fig. 13 with SU{N) boson exchange vanishes, because TrT"^ = 0. The 
net result is that the difference between SU{N) and U{N) is order 



Problem 3.1 

Prove eqs. (3.4), (3.5). 
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Fig. 14. The quark propagator in terms of the one-particle irreducible self-energy S. 



3.2. The 't Hooft Model 

The 't Hooft model is large N QCD in 1 -|- 1 dimensions. This theory 

was solved by 't Hooft [1,2] to obtain the exact meson spectrum. There is 
an extensive discussion of this model in Coleman's lectures [3]. I will not 
repeat the solution of the model here. Instead, I will summarize why the 
model is solvable, and show some of the numerical solutions. 

In 1 + 1 dimensions, the gauge coupling constant g has dimensions of 
a mass, and becomes relevant in the infrared. The theory is confining, 
and has a linear potential at large distances. Even QED confines in 1 + 1 
dimensions, and has a linear potential. It is convenient to use light-cone 
coordinates 

with metric 

ds^ = 2dx~^dx~ , g^v = 

and light-cone gauge ^+ = A- = 0. The field-strength tensor has a single 
non- vanishing component, 

F+_ = d+A_ - d-A++i-^ \A+,A_] = -d-A+, 

so the theory becomes effectively Abelian. This is the first important sim- 
plification. The second is planarity, which allows one to solve the theory 

exactly. 

Consider the quark propagator, which can be written in terms of the 
one-particle irreducible piece S, as in fig. 14. The equation for S is given 
graphically in fig. 15. Planarity is crucial for the derivation of this relation. 
We earlier derived the iV-counting rules only for vacuum graphs. The 
results can easily be extended to the fermion two-point function, which is 
obtained by differentiating a vacuum graph once with respect to a fermion 
bilinear source, i.e., by cutting the closed quark loop at one point. Planarity 
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Fig. 15. The equation for S. 





+ 



+ 



+ ... 



Fig. 16. Rainbow sum for the exact quark propagator. 



of the vacuum graphs then impHes that the quark propagator is given by 
planar graphs with all gluons on one side of the qnark line. The first gluon 
emitted by the quark must be the last gluon absorbed, because the diagram 
is planar, and there are no gluon self-interactions. This immediately leads 
to fig. f 5. 

The equation for S in fig. 15 can be iterated, to give the solution fig. 16 
for the quark propagator. The diagrams in fig. 16 are known as rainbow 
diagrams, and the rainbow approximation is exact in the 't Hooft model. 
The analytical solution of fig. 15 is straightforward [3]. The solution is that 
the quark propagator is given by the free-quark form 1/ (|< — A/) , with the 
renormalized quark mass M given in terms of the bare quark mass by 



Note that there is still a pole in the quark propagator, even though the 
theory is confining, and there are no quark states in the spectrum of the 
theory. Also note that the pole in the quark propagator can be tachyonic 
if < (/^/27r. Nevertheless, the theory has a sensible meson spectrum for 
all values of > 0. 

The meson propagator has the exact graphical expansion in fig. 17, and 
the ladder graph approximation is exact in the large N limit. This again 
follows trivially from the planar diagram structure of the theory and the 
absence of gluon self interactions: gluons are not allowed to cross in any 
graph. The Bethe-Salpeter equation for the meson wavefunction, known as 
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+ 



Fig. 17. Graphical expansion for the meson propagator. The shaded blobs are full quark 
propagators. 



Fig. 18. The Bethe-Salpeter equation for the meson wavefunction. The cross-hatched 
blob is the meson wavefunction <i>(x). 



the 't Hooft equation, is shown graphically in fig. 18, and follows from the 
meson propagator fig. 17. Let P be the total momentum of the meson, q be 
the momentum of the quark (the quark is well-defined, since there is no pair 
creation in the large N limit), x = q^/P^, and 4>{x) the amplitude to find 
the quark with this light-cone momentum fraction. A simple calculation 
leads to the 't Hooft equation [1-3] for the meson wavefunction 



where P denotes the principal value, /U is the meson mass. Mi is the renor- 
malized quark mass, and M2 is the renormalized antiquark mass. 

This equation can be solved numerically using a Multhopp trans- 
form [10]. The ground state wavefunction (j){x) for m\ = 1712 = 1 (renor- 
malized masses in units of g/\/27r) is shown in fig. 19. The meson mass 
is n = 2.7, which is larger than the sum of the two quark masses. The 
wavefunction of the first excited state is shown in fig. 20, and corresponds 
to a meson with mass fj, = 4.16. One can show that for large excitation 
number n, the meson mass is linear in n. The ground state wavefunction 
for mi = 1712 = 10 is shown in fig. 21; the meson mass is /i = 20.55. Clearly, 
increasing the quark mass narrows the momentum spread in the wavefunc- 
tion, and also decreases — mi — m2, the strong interaction contribution 
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0,3 0,4 0,6 0,8 



X 

Fig. 19. The ground state wavefunction of the 't Hooft model with mi = m2 = 1. 
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0,3 0,4 0,6 0,8 

X 



Fig. 20. The wavefunction for the first excited state of the 't Hooft model with mi = 

7712 = 1. 



to the mass. For unequal masses, the momentum distribution is asymmet- 
ric, and the heavy quark carries most of the momentum. Figure 22 shows 
the ground state wavefunction for mi = 5, m2 = 1, with a meson mass 
H = 6.72. For hght quarks, m ^ 0, the meson wavefunction is not affected 
much by the quark mass. As one might expect, the structure of the wave- 
function is governed by the scale g rather than m. The wavefunctions of 
the 't Hooft model show many of the features one would expect for the 
wavefunctions of mesons in QCD in 3 -|- 1 dimensions. 
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0,2 0,4 0,6 0,8 1 



X 

Fig. 21. The ground state wavefunction of the 't Hooft model with mi = m2 = 10. 




0,2 0,4 0,6 0,8 



X 

Fig. 22. The ground state wavefunction of the 't Hooft model with mi = 5 and m2 = 1. 



3.3. N -Counting Rules for Correlation Functions 

We have discussed the 7V-eounting rules for connected vacuum diagrams. 
These can be used to derive A^-counting rules for gluon and quark cor- 
relation functions. The correlators we will study are vacuum expectation 
values of products of gauge invariant quark and gluon operators. The op- 
erators need not be local; all we require is that they cannot be split into 
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Fig. 23. Leading contribution to a quark correlation function. The shaded region repre- 
sents planar gluons, and (S) are insertions of the fermion bilinears Oj. 



pieces which are separately color singlets. Allowed operators are 
but not 

Operators involving quark fields must be bilinear in the quarks. Let Oi 
denote a generic operator written in terms of rescaled fields, and add the 
source term NJiOi to the rescaled Lagrangian eq. (3.7). The entire La- 
grangian still has an overall factor of TV, so the A^-counting rule eq. (3.8) 
still holds. Connected correlation functions are then obtained by differen- 
tiating W{J), the sum of connected vacuum graphs, with respect to the 
sources 

The order N'^ contribution to W{J) is from graphs with only gluon lines. 
This can contribute to the correlation function {O1O2 ■ ■ ■Or)c^ provided 
none of the operators contain any quark fields. Thus pure-gluon r-point 
correlation functions arc of order N'^~^ . The diagrams that contribute to 
these are planar diagrams with insertions of Oi. The first contribution to 
W{J) that involves quarks is of order N. Thus r-point correlation functions 
that involve quark fields arc of order N^~'^ . The diagrams that contribute 
to them are planar diagrams with a single quark loop with the fermion 
bilinears inserted on the quark line, as shown in fig. (23). 

The A'^-counting rules were derived using the rescaled Lagrangian 
eq. (3.7). One can obtain A^-counting rules for correlation functions with 
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fields normalized as in the original Lagrangian eq. (3.1) by using the re- 
placement 



ip = \/Nip, A 




it particular, quark bilinears are related by a factor of A'', Oj = NOi, e.g. 
ipijj = Nijjip. 

The A'^-counting rules for correlation functions can be used to derive 
the TV-counting rules for meson and glueball scattering amplitudes. We 
will use the notation Gi to denote a gauge invariant pure-gluonic operator, 
and Hi to denote a gauge invariant operator bilinear in the quark fields, 
where the operators are written in terms of the rescaled fields. Gluon op- 
erators Gi can create glueballs, and quark bilinears Hi can create mesons. 
The two point function {GiG2)c is of order imity, so Gi creates glueball 
states with unit amplitude. The r-point function (Gi . . .Gr)c is of order 
N'^~'^. Thus r-glueball interaction vertices are of order N'^~'^, and each 
additional glueball gives a 1/N suppression. The meson two point corre- 
lation function {HiH2)c ~ l/N. Thus VNHi creates a meson with unit 
amplitude. The r-point function {\/NHi . . .\/NHr)c is of order N'^~'^/^. 
Thus r-meson interaction vertices are of order N'^~^/'^, and each additional 
meson gives a 1/VN suppression. Finally, one can look at mixed glue- 
ball meson correlators, (Gi . . . Gr-ffi . . . Hs)c, which are of order N^~'^~^. 
Thus an interaction vertex involving r glueballs and s mesons is of order 
so that (GVnH) is of order 1/Vn — meson-glueball mixing is 
of order l/ViV, and vanishes in the N oo limit. 

One important result that we will need later is that the pion decay 
constant /jr is of order V^. The matrix element (0|g7''75T°g|7r''(p)) = 
ifTrP^S"^, where the quark fields are normalized as in the original La- 
grangian eq. (3.1). The A^-dependence of the matrix element can be ob- 
tained from {Hi^/NH2) ^ 0{1/Vn), where the first bilinear is the axial 
ciirrent, and the second bilinear produces a pion from the vacuum with am- 
plitude of order unity. One needs to multiply Hi by an additional factor 
of N to convert the axial current from rescaled fields to the original quark 
fields, so that 

/.oc ViV. 

The iV-counting rules imply that one has a weakly interacting theory of 
mesons and glueballs with a coupling constant As in any weakly 

interacting theory, one can expand in the coupling constant l/\fN. The 
leading order graphs are tree-graphs, and the leading order singularities 
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are poles. At one-loop (i.e. one gets 2 particle cuts, at two loops, 

three-particle cuts, and so on. QCD, a strongly interacting theory of quarks 
and gluons, has been rewritten as a weakly interacting theory of hadrons. 
The leading (in N) interactions bind the quarks and gluons into color sin- 
glet hadrons. The residual interactions between these hadrons are 1/N 
suppressed. The 1/N expansion is also equivalent to the semiclassical ex- 
pansion for the meson theory. These results will also hold for baryons. 

The spectrum of the theory contains an infinite number of narrow glue- 
ball and meson resonances. The resonances are narrow, because their decay 
widths vanish as iV ^ oo, since all decay vertices arc proportional to pow- 
ers of the weak coupling constant and hadron masses (i.e. phase 
space factors) do not grow with N. There must be an infinite number of 
resonances to reproduce the logarithmic running of QCD correlation func- 
tions. A meson two-point correlation function can be written as a sum over 
resonances, 

/d4xe^«-(g(x)g(o))^ = ^-^^, 
J i 1 

since single meson exchange dominates in the large N limit. The left hand 
side has logarithms of q^, which can only be reproduced by the right hand 
side if there are an infinite number of terms in the sum. 

3.4. The Master Field 

The A^-co\inting rules imply that the functional integral measure is concen- 
trated on a single gauge orbit, and that fluctuations vanish in the N ^ co 
limit. For example, if G is a gauge invariant operator made of gluon fields, 
{G) 0{N), and its variance is 

Thus AG/G ~ 0{1/N) 0. It is easy to see that all gauge invariant 
observables have no fluctuations in the N ^ 00 limit. The functional inte- 
gral measure must then be concentrated on a single gauge orbit, known as 
the master orbit, represented by a set of gauge equivalent vector potentials 
^^(a;). It is expected that one can find a point on this orbit where the 
vector potential is given by constant matrices A^. This is the master field, 
and all correlation functions are given by computing them in this single 
field configuration. A lot of information can be encoded in the master 
field, since it is an 00 x 00 matrix. The master field has recently been 
determined for QCD in 1 -|- 1 dimensions [11,12]. 
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CO 




Ci 



C2 



Fig. 24. Examples of Wilson loops. Ai and A2 are the unshaded and shaded areas, 
respectively. 



There are some other interesting results for large N QCD in 1 + 1 di- 
mensions which have been obtained recently. The basic gauge invariant 
observables in a gauge theory are the Wilson loops. Kazakov and Kostov 
obtained the exact expression for all Wilson loops [13]. The Wilson loop 
for a close curve C is the expectation value of the trace of the path-ordered 
exponential 



and is of order N. For a simple closed curve, the Wilson loop for 1 + 1 
dimensional large N QCD satisfies an exact area law, 



where A is the area enclosed by C. More interesting examples are self- 
intersecting curves, such as those in fig. 24, with Wilson loops 



W{C\) = TVe-s'^-^^+^^'Z^^ W{C^) = 7Ve-9'(^^+2^=)/2 [l - , 



where Ax and A-i, are the unshaded and shaded areas, respectively. The 
master field for 2D QCD has been shown to reproduce the Kazakov-Kostov 
results for the Wilson loops. 

Problem 3.2 

In a non-Ahclian theory, the field strength tensor does not determine all 
the gauge invariants. In 1 + 1 dimensions for gauge group SU (2), construct 
two different vector potentials which both produce the same constant field 
strength tensor F12 = fr^, where f ^ is a constant, and yet give different 
values for Wilson loops. 




(3.11) 



W{C)=Ne-s"^/^, 
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Problem 3.3 

DeGnc the matrix 

Ur = Pe-'^ ■1'^'^ 

where = Af^T"^ with in representation R of SU{N). The Wilson 
loop in the fundamental representation (denoted by F), eq. (3.11), is 



where (•) represents a functional integral over all gauge held configurations. 
Show that the Wilson loop in the adjoint representation is 



This equation holds in any number of dimensions, and does not depend on 
taking the large N limit. The Wilson loop W>(C) is expected to have an 
area law, 

Wf{C) ~ N exp(-A area(C)), 

where A is an unknown constant. Show that in the large N limit, one 
expects [14] 

WadjiC) ^ cxp(-2 A arca(C)) +cxp(-A' pcrimctcr(C)), 

where A' is another unknown constant. What are the implications of this 
formula for confinement /screening of adjoint quarks? 

4. Meson Phenomenology 

4.1. Zweig's Rule 

Meson correlation functions at leading order in 1 /N are given by diagrams 
with a single quark loop. Annihilation graphs, such as those in fig. 25, have 
two quark loops and arc suppressed by one power of A'^. The suppression 
of annihilation graphs is known as Zweig's rule. One consequence of this is 
that mesons occur in nonets for three light quark flavors. There are nine 
possible qq states, which are divided into flavor octets qT^ q, where T^ is a 
SU (3) flavor matrix, and a singlet qq. Usually, the singlet and octet mesons 
are not related, because the singlet meson can mix with gluons. In the large 
N limit, this mixing is suppressed. The singlet and octet meson couplings 
are related by treating them as members of a U{N) multiplet, with the 



Wf{C) = {TiUf) 
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Fig. 25. An annihilation graph, which violates Zweig's rule and is l/N suppressed. 



Fig. 26. The leading order graphs for the ^(gg)^ (<l<l)i correlation function. The shaded 
regions represent planar gluons. 



[/(I) and SU{N) couplings having the same normalization. Examples of 
this are: in the large N limit, the vector meson nonet of {p,u},(l),K*} all 
have the same mass, and the decay constants of the tt and rj' are equal, 



4-2. Exotics 

At leading order in l/N, there are no exotic states. This does not mean that 
there arc no exotic states in QCD, only that their binding is sublcading in 
the l/N expansion. As an example, consider the propagator {{qq}^ {qq)^} 
for a four-quark state qqqq. It is easy to see that at leading order, the 
graphs which contribute have the form shown in fig. 26, i.e. the correlation 
function is {qq qq)^ , the square of the qq correlation function. This is the 
correlation function for two non-interacting mesons, and so {qq)^ creates 
two mesons, rather than a four-quark state. 

4-3. Chiral Perturbation Theory 




Jtt — fri'- 



The U{Nf)l X U{Nf)r chiral symmetry symmetry of QCD is spon- 
taneously broken to a diagonal U{Np)v vector symmetry, resulting in 
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(pseudo-) Goldstone bosons.* This form for the breaking can be proved in 
the large hmit [15]. The low-energy interactions of the pseudo- Goldstone 
bosons of QCD, the tt, K, rj, and r]', can be described in terms of an effec- 
tive Lagrangian known as the chiral Lagrangian. One can imagine comput- 
ing the chiral Lagrangian by evaluating the QCD functional integral with 
sources for the pseudo-Goldstone bosons. The source terms are fermion bi- 
lincars. In the large N limit, we have seen that the leading order diagrams 
that contribute to correlation functions of fermion bilinears are of order 
N, and contain a single quark loop, as in fig. 23. This implies that the 
leading order terms in the chiral Lagrangian arc of order N. It is also clear 
from the structure of fig. 23 that the leading order terms can be written 
as a single fiavor trace, since the outgoing quark flavor at one vertex is 
the incoming flavor at the next vertex. Similarly, diagrams with two quark 
loops have two flavor traces, and are of order unity, and in general, those 
with r quark loops have r traces, and are of order N^~'^. 

The chiral Lagrangian is written in terms of a unitary matrix 

U = e'^'^/f\ (4.1) 

where /■„ w 93 MeV is the pion decay constant, and 



\/2 



(4.2) 



is the matrix of pseudo-Goldstone bosons. The r]' has been included, since 

it is related to the octet pseudo-Goldstone bosons in the large N limit, by 
Zweig's rule. The order terms in the chiral Lagrangian are 

i(2) = ^^Tr D^'UD^U-'^ + ^BTr {m^U + mU-'^) , (4.3) 

where m is the quark mass matrix in the QCD Lagrangian. The first 
term is order N, since /^r oc ^/N. The second term in eq. (4.3) also has 
a single trace and is of order A^, so -B is of order unity. The U field has 
an expansion in powers of tt/Ztt. Thus each additional meson field has a 
factor of I/./tt oc 1/Vn, which gives the required suppression for 

mesons derived earlier. The effective Lagrangian eq. (4.3) has an overall 
factor of N, and the U matrix is independent, so the 1/N expansion is 

The axial anomaly is O (l/N). See section 4.6. 
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equivalent to a scmiclassical expansion. Graphs computed using the chiral 
Lagrangian have a suppression for each loop. 

The order terms in the chiral Lagrangian are conventionally written 
as [16] 

L^^^ =Li \TvD^U-'^D'^uf + L2Ti:Df,U-'^D^U'Tr:D^'U-'^D''U 

+L3 Tr Df,U-^Di'UD^U-^D''U 

+L4. Tr D^U-^Di'U Tr {U-^m + w)U) 

+L5TvD^U-^DfU {U-^m + m)U)+LQ [Tr {U-^m + m)U)]^ 
+Lr [Tr {U-^m - m^U)] % Tr {m^Um^U + U-'^mU-'^m) 

+L,oTrU^F^''U-^FL^,, (4.4) 

where Fl,r are the field-strength tensors of the (flavor) U{3)l and J7(3)ij 
gauge fields. The terms with a single trace, L3, L5, Lg, Lg and Liq should 
be of order N, and those with two traces, Li, L2, L^, Lq and L7 should 
be of order one. This is not correct, because of one subtlety. There is an 
identity 

Tr ABAB = -2 Tr A^B^ + 1 Tr Tr + (Tr ABf , (4.5) 

which is valid for arbitrary traceless 3x3 matrices A and B. Using A = 
Dfj,UU~^ and B = D^UU^^ in eq. (4.5) gives the relation 

Tr D^UDM~'^Di^UD''U-^ = -2Tv D^UD^'U-^D^UD^'U-^ 
+ - Tr DMD^'U-'^ Tr DMD^U-'^ +Tr DJID^U''^ Tr D^UD^U''^. 



It is important to remember that this relation is special to three flavors, 
and would not hold for an arbitrary number of flavors. The operator 
Tr D^UDvU~^D^UD'^U~^ is a single trace operator and can occur in the 
Lagrangian with a coefScient c which is of order N. Eliminating the oper- 
ator using the identity eq. (4.6) gives the contributions SLi = c/2, 6L2 = c 
and 6Ls = —2c to ii_3. L3 was already of order A'', and so remains order 
A''. Li and L2 arc now of order N, because of the c term, but 2Li — L2 is 
of order unity. Thus one finds the A^-counting rules 



2 



(4.6) 



0{N) : 
0{1): 



Li, L2, L3, L5, Ls, Lg, Lio 
2Li — L2, L4, Lq, L7 
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Table 1 

Experimental values for the coeffi- 
cients of the order terms in the 
chiral Lagrangian, eq. (4.4). Values 
are from ref. [18] 



Li 


Value 


Order 


2Li - L2 


-0.6 ±0.5 


1 


u 


-0.3 ±0.5 


1 


Le 


-0.2 ±0.3 


1 


Lr 


-0.4 ±0.2 


1 


L2 


1.4 ±0.3 


N 


Lz 


-3.5 ± 1.1 


N 


L5 


1.4±0.5 


N 


Ls 


0.9 ±0.3 


N 


Lg 


6.9 ±0.7 


N 


Lio 


-5.5 ±0.7 


N 



These are the A^-counting rules given in ref. [16], with the exception of 
L7, which is taken from ref. [17]. In ref. [16], Lj was argued to be of 
order A''^. We will return to Lr in section 4.7, after discussing the 77'. The 
experimental values for the L's are given in Table 1. The terms of order N 
are systematically larger than those of order unity. 

Higher derivative terms in the chiral Lagrangian are suppressed by pow- 
ers of the chiral symmetry breaking scale ~ 1 GeV [19,20]. In the large 
A'' limit, is of order unity, and so stays at around 1 GeV. Loop graphs in 
the chiral Lagrangian are proportional to l/liirf^)^ and are of order 1/A^. 
Thus in the large N limit, the chiral Lagrangian can be used at tree level, 
and loop effects are suppressed by powers of 1/N. 

4-4- Non-leptonic K Decay 

Weak decays of hadrons can be computed in terms of an effective low 
energy Lagrangian, since hadron masses are much smaller than the W 
mass. Semileptonic weak decays, such as K ^ ttIv, can be computed 
using the weak Lagrangian 

where = (1 — 75)72 is the left-handed projection operator. To all orders 
in the strong interactions (and neglecting electromagnetic corrections), the 
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(a) (b) (c) 

Fig. 27. Diagrams for K — » tttt. There are related diagrams with K ^ w. The shaded 
regions represent planar gluons. The wavy line is the W boson. The (S) are insertions of 
the axial current. Diagram (a) is of order N'^. The factorization violation contributions 
(b) and (c) are order N and order unity, respectively. 



matrix element for the decay Hi — > HfiU can be written as 

-i^ {Hflu^^PLS \Hi) {eul h^PLV |0) , 

i.e. it factorizcs into the product of a hadronic matrix clement, and a Icp- 
tonic matrix element. The leptonic part can be computed explicitly using 
free Dirac spinors. The hadronic part is the matrix element of a current, 
and for K decays, it can be computed in terms of the decay constant /k- 
The decay constant is determined from the measured K fiU decay rate, 
and can then be used to predict other decay rates such as ii' — > -KeV, etc. 

Non-lcptonic K decay amplitudes are more difficult to compute, since 
they depend on the hadronic matrix elements of four-quark operators. To 
analyze these using the large A'' limit, it is convenient to first look at the 
weak amplitudes directly in terms of W exchange, rather than using the 
effective weak Lagrangian. The K ^ tttt amplitude in the large N limit, 
and to lowest order in the electroweak interactions, is given by diagrams 
with a single W boson. The W boson is color neutral, so the A^-counting 
for a diagram with a W boson is the same as that for the diagram with 
the W boson removed. The diagram must contain a quark loop, since the 
K and tt mesons contain quarks. The leading order diagrams are planar 
diagrams with two quark loops (fig. 27(a)). These are disconnected dia- 
grams when the W boson is removed, so each quark loop subgraph is of 
order N, and the overall graph is of order N'^. The K ^ tttt amplitude 
is order -x/iV, since each current produces a meson with amplitude ^N. 
This should be compared with a three-meson amplitude in the strong in- 
teractions, which is l/\fN. Weak interaction perturbation theory is an 
expansion in Gf^qctiN . Formally, this diverges if one takes N ^ oo with 
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Gf fixed. This docs not mean that one should add lots of W's to strong 
interaction processes to get an amplitude that grows with N. In the end, 
the results are going to be applied to A'' = 3, with Gp set to its experimen- 
tal value. One can first use perturbation theory in Gi? to write the weak 
decay amplitudes in terms of hadronic matrix elements, and then apply the 
1/N expansion to evaluate the purely strong interaction matrix elements. 

There are no gluon exchanges between the two quark loops, so the lead- 
ing order K ^ mr amplitude fig. 27(a) has clearly factorized into the 
product of (0| \K) x (ttttI j^vy |0) plus terms with -k ^ K, where is 
the weak current. Factorization is exact in the large A'^ limit. Corrections 
to the factorization approximation, such as figs. 27(b,c) are suppressed by 
1/A'' and l/N"^, respectively. The K ^ tttt amplitudes can be computed 
in terms of fx in the factorization approximation, since each hadronic ma- 
trix element is that of a current. In particular, the ratio ^1/2/^3/2 of the 
A7 = 1/2 and A7 = 3/2 amplitudes is equal to \/2. One easy way to 
compute this is to note that the amplitude for tt'^tt^ from fig. 27(a) 

vanishes, since all the mesons are neutral, and the W boson is charged. 

Experimentally, ^1/2/^3/2 ~ 21, which is the famous A/ = 1/2 rule in 
non-leptonic K decays. There is no A/ = 1/2 enhancement at leading order 
in the 1/N expansion [21,22]. At first sight, this is a disaster. However, it is 
important to keep in mind that non-leptonic weak decays are a multiscale 
problem, and involve both M\y and Aqcd- Renormalization group scaling 
of the effective weak Lagrangian from Mw to a low scale /u ~ 1 GeV 
produces an enhancement of the A7 = 1/2 amplitude. Formally, this 
enhancement is 1/N x log Mw /m- = 3 and log Mw //i ^ 4, so one should 
sum all powers 1/Nx log Mw / This is done by using the renormalization 
group to scale the weak Hamiltonian down to some low scale n of order the 
the hadronic scale. Matrix elements of the low-energy weak Hamiltonian 
do not contain any large logarithms, and should have 1/A^ corrections of 
canonical size. This has been examined in detail, and is claimed to produce 
a satisfactory understanding of the A/ ~ 1/2 nde [23]. The analysis is 
involved and will not be repeated here. A simpler example is considered in 
the next section. 

Problem 4.1 

The K decay amplitudes are 



A{K+ 




A {K° TT+TT 



) 



Al/2 + -1^^3/2 
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(a) 



K 



(b) 



Fig. 28. K — K mixing diagrams. The blob represents the four-quark operator 
d'yf^ Pls d'y'^ Pl^s, and the shaded regions represent planar gluons. The operator ver- 
tex has been made "transparent" so that the color flow along the fermion lines is visible. 
Diagram (a) is of order N'^, and Diagram (b), is of order N. 



neglecting final state interaction phases. Compute the amplitudes 
TT+TT^, K'^ 7r"'"7r~, and K'^ 7r°7r° in the factorization approximation, 
and use these to obtain and Ai/2- Compute the decay widths for 
and Kg decay, and compare with experiment to obtain and Ai^2- 
Note that K%^K°. 

4.5. K — K mixing 

The K — K mixing amplitude is of second order in the weak interactions. 
The mixing amplitude is given by the matrix element of the AS = 2 La- 
grangian between K and K. The AS* = 2 Lagrangian is 

L = ^riifj) d^^^PLsd^PLS, (4.7) 

where 77 has dimension two, can be computed in perturbation theory, and 
includes renormalization group scaling from Myy down to some hadronic 
scale /U. The ^ dependence of r] is cancelled by the anomalous dimension 
of the four-quark operator. It is conventional to write the K — K matrix 
element of the four-quark operator as 

{K\ dji'PLsdrPLS \K) = yiMlBK{li), 

where BKifj) parameterizes the hadronic matrix element of the four-quark 

operator renormalized at ^. In the large A'' limit, the matrix element is 
given by fig. 28(a), and factorizes into the matrix element of two currents, 

(K\ d-yf^PLsd^PLS \K) = 2 (K\ d'y'^PLS |0) x (0| I-j'^Pls \K) . 
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The factor of two arises because there arc two ways of combining the weak 
currents with the mesons. Corrections to the factorization approximation 
are of order 1/N from fig. 28(b). The matrix element (0| d'y^^P^s = 
IkP^/V^^ so in the factorization approximation, one finds Bk = 3/4 + 
0{l/N) [24]. The same result holds for B-B mixing, = 3/4+e> (l/A^). 
There is no ^ dependence in B in the large A'' limit, because the anomalous 
dimension of the four-quark operator is of order 1/N . The 1/N corrections 
to Bk,b do not contain the scale Mw, since renormalization group scaling 
has already been used to obtain the effective interaction eq. (4.7) at a low 
scale. One does not expect enhanced 1/N corrections in Bk.b- Recent 
lattice results give Bk{2 GeV) = 0.62 ±0.02 ±0.02 [25] and Bk{2 GeV) = 
0.628 ± 0.042 [26], which are consistent with 3/4 + O {1/N). 

4-6. Axial U{1) and the r/ Mass 

The U{1)a flavor symmetry is broken by anomalies, so the rj' is not a Gold- 
stone boson. The anomaly graph involves a quark loop, and is suppressed 
in the large N limit. This allows us to study the r]' as a pseudo-Goldstone 
boson, with 1/N as a symmetry breaking parameter. 
The QCD Lagrangian including the 9 term is 

L = -i Trf^.i^"- + £2^^ TrF^^F^'^+iiiip - m) ^, (4.8) 

where the usual coupling constant has been replaced by /N . is a pe- 
riodic variable with periodicity 27r. It is clear from the form of eq. (4.8) that 
quantities will depend on the combination 6/N, since that is the parameter 
combination which occurs in the QCD Lagrangian. 

There is no ^-dependence of any physical quantity in perturbation the- 
ory, which makes the analysis of the U{1)a sector tricky. The vacuum 
energy E is of order iV^ . If one assumes that there is no 1 /N suppression 
of the 9 dependence, it must have the form 

E = N^F{9/N), 

where F is some function with periodicity 2Tr/N. In the dilute instanton 
gas approximation, the 9 dependence of quantities has the form 

in the one-instanton sector. This is exponentially suppressed in N, and one 
might think that all 6 dependence is exponentially small in N. This con- 
clusion is believed to be incorrect. The dilute instanton gas approximation 
is not valid because of infrared divergences. 
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There is one important result about the 9 dependence of QCD when 
fermions are included — if any fermion is massless, all 9 dependence van- 
ishes. This leads to the following puzzle: fermion loop contributions to the 
vacuum energy are order N, so how can they cancel the order N'^ vacuum 
energy of the pure glue theory? To study this, look at the second derivative 
of the vacuum energy with respect to 9 [27,4] , 



d9^ ~ [ 8n'^N 



j d'^x (q T {tv FF{x)T:v FF{<S)^ Oj, (4.9) 



which is 0{1). Define 

U{k) = j d^x e'*^-^ (o T (tiFF{x) TiFF{0) 



so that 

fE 
d9^ 



9 



U{0). 



One can write the two-point correlation function as a sum over intermediate 
single-particle states. 



U{k)= J2 



glueballs 



E 



Nbi 



where we have used the A/'-counting rules (0|Trf Fjglueball) N, 
(0|Tri^F|meson) ~ V^. Multiparticle states can be neglected in the large 
A'' limit. Here to„ and M„ are glueball and meson masses, and Nun and 
n are the amplitudes for TrFF to create a glueball or meson from 
the vacuum. In the pure- glue theory, only the first term is present, so 
U{0) ~ N'^ and d^E/d9'^ ~ 1. In the theory with quarks and gluons, the 
second term is also present. The only way that the second term can cancel 
the first is if one meson has a mass of order 1 / \fN . This can cancel the 
first term at fc = 0, which is what is needed to cancel the 9 dependence of 
the vacuum energy, but does not cancel the first term at arbitrary values 
of k. It is believed that the rj' mass is of order l/\/]V, and produces the 
required cancellation.* With this assumption, one finds that 



no quarks 



N 



One might wonder how two terms of the same sign can cancel ea<;h other, 
lution is that there is an equal time commutator that must be added to eq. 
the appendix of ref. [4] for more details. 



The reso- 
(4.9). See 
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The matrix element VNbr,'g'^ /Sn'^ = {g'^ /8Tr'^){0\Tr FF\r]'), which can be 
written as 

(^VStt^) (o| TtFF \rj') = ^ (0| |r?') = ^^M^^, 
using the anomaly equation 

This gives the Veneziano-Witten formula for the rj' mass, 

2 27Vf / d^i^no quarks \ . . 

" U V ^^^^ ^=0 ^ ' 

The 77' is a Goldstone boson in the large A'' limit, and M^, is linear in the 
symmetry breaking parameter 1/N. In general, one can show that the rj' 
dependence of a zero-momentum amplitude in the theory with quarks can 
be obtained from the 9 dependence in the theory without quarks, by the 

replacement 

9^9+^2N^r^'/U. (4.11) 

The f]' mass is (PE/drj'"^ = M^y, which reduces to eq. (4.10) using cq. (4.11). 

The form eq. (4.11) can also be obtained using the U{1)a Ward identity. 
Under an axial U(l) transformation, tpL e^'^ipL, V'-R ^ e~'"f/'-R) ^^'i 
9^9 — 2Npa. The 77' is the Goldstone boson of the U(1)a symmetry, 
and transforms as rj' ^ rj' + f^ra, so the U{1)a invariant combination is 
eq. (4.11). 

In chiral perturbation theory, the U matrix can be extended to include 

the 77', 

where = at leading order in l/N* Then the linear combination in 
eq. (4.11) is 

6* - i log det U. 

One can obtain the zero-momentum rj' couplings by using this linear com- 
bination for all the 9 dependence in the effective Lagrangian. There 

arc also momentum-dependent rj' interactions which are not related to ^- 
depcndence. and cannot be determined by this method. 

This has already been done in the U of eq. (4.2). 
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In the largo N and chiral limits, tlio tt. K, rj, r]' nonet is masslcss. Includ- 
ing non-zero quark masses and corrections gives mass to the mesons. 
For simplicity, consider m„ = md = 0, rUg ^ 0. The neutral mesons in 
the nonet can be chosen to be mu, dd and ss, rather than the tt*^, t] and 
rl . Since m„ = rrid = 0, there is an exact U(2) x ?7(2) chiral symmetry in 
the large N limit, so the meson mass matrix for the uu, dd and ss mesons 
must have the form 



= , (4.12) 



where C is some function of m^. There can be no off-diagonal terms, since 
Zwcig's rule is exact in the large N limit. The 1/N correction to the mass 
matrix has the form 





(4.13) 



(in the limit of equal quark masses), since the amplitude for qiqj^ If'^lf 
does not depend on the flavors i, f of the initial and final quarks. The 1/A'' 
factor is explicit, so that a, which represents the strength of the annihilation 
graphs, is of order one. The quark mass and 1/N are both treated as small, 
so that effects of order niq/N have been neglected. The complete mass 
matrix is the sum of eqs. (4.12)-|-(4.13). It has one zero eigenvalue (the 7r°) 
since chiral SU{2) x SU{2) is still an exact symmetry. The two non-zero 
eigenvalues give the ratio [28] 



M2 3 + i? - V9-'2R + i?2 



R = CN/a. (4.14) 



M2 3 + i? + \/9 - 2i? + i?2 ' 
Irrespective of the value of R, one finds 

The experimental value for the ratio is 0.33, which exceeds the bound, but 
not by much. (Remember that our expansion parameter is about 1/3.) 
The bound was derived neglecting m,j_d, and keeping the lowest order term 
in 1/N. Including light quark masses, and adding the 1/A^ correction gives 
a bound that is consistent with the experimental value [29] . 

The rj' mass M^^, is a function of the quark masses and 1/N. In the 
limit — » 00, it is of order uig, and in the chiral limit niq 0, it is of 
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order so it is sometimes said that the large N and chiral Hmits do not 
commute. The origin of this non-commutativity is clear; the mass matrix 
is the sum of two terms eqs. (4.12)+(4.13), and the eigenvalues depend on 
the ratio z ~ Nrriq/A-^, where is a typical hadronic scale which has been 
used to make z dimensionless. The large N limit is z ^ oo, and the chiral 
limit is 2; — > 0, and M^, has different forms in these two limits. However, 
what is relevant for applying large; and chiral symmetry is that 1 /N and 
niq/A^ are both small; their relative size z is irrelevant. 

Chiral perturbation theory is an expansion in derivatives over some chi- 
ral symmetry breaking scale A^, which is 0(1) in the large N limit. The 
rj' is light in the large A'^ and chiral limits, irrespective of the relative size 
of 1/7V and niq. Since it is light, the rj' should be included as an explicit 
degree of freedom in the large N chiral Lagrangian, to avoid an inconsistent 
expansion. The U{Nf) x U{Nf) chiral Lagrangian in the large N limit to 
order has been worked out in ref. [30] 

Problem 4.2 (The rj' Mass) 
Derive eqs. (4.14) and (4.15). 

4-7. Resonances and 1/N 

In section 4.3, we derived the N dependence of an effective Lagrangian for 
mesons. The leading order terms in the Lagrangian are of order N, and 
have a single trace over flavor. The first correction is of order unity, and has 
two flavor traces, etc. In addition, every meson field carries a suppression 
factor of 1 /VN. The Lagrangian can be represented schematically as 

L = NTrX + TrYiTrY2 + ^TrZiTrZ2TrZ3 + (4.16) 

where X . Yi and Z,; are functions of M /\fN , where M is a meson field. Here 
Tr Yi Tr Y2 is an abbreviation for the sum of all possible terms written as 
the product of two traces, etc. For example, in the chiral Lagrangian, Tr AT 
represents terms such as D^U^^D'^U , where U = exp(2i7r//), with / of 
order Here we will consider a more general low energy Lagrangian, 

that includes the Goldstone bosons as well as additional meson fields, and 
study the form of terms induced by integrating out heavy meson fields. 

In the large N limit, mesons form nonets. It is convenient to represent 
a meson nonet (such as the p, lj, (f) and K*) by a 3 x 3 matrix M. The 
one-meson couplings can be obtained from the Lagrangian cq. (4.16), by 
retaining the terms which contain one power of Mj^pN , and schematically 
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have the form 



= ^/NTrAM + ^TrBMTrC 



+ 



N 



TTDMT1E1T1E2 + .. 



(4.17) 



The terms induced by integrating out the meson multiplct M at tree level 
can be obtained from eq. (4.17). The meson propagator (at zero momen- 
tum) is 



1 



o9 r69 



mi ' 



-^S''^ a, 6=1,.... 

-KS"'' 0,6 = 9 



(4.18) 



where mg is the mass of the octet mesons, 

Ol/m2 = 

is related to the mass difference of the the octet and singlet mesons, and 
the first S function in eq. (4.18) is over a, 6 = 1, . . . , 9. Writing M as M°'T°- 
(a = 1-9), and using the identity eq. (3.5) with A'' replaced by Np = 3, 
shows that the terms induced by meson exchange at lowest order in the 
derivative expansion are 

\/7V Tr ^T" + Tr BT" Tr C + Tr DT" Tr £1 Tr £2 



X A"*" X 



TV Tr AT 
1 



1 



Tr BT" Tr C 



1 



2mi 



1, 

1 

N 



NTrAA + 2TiABTrC + 
1 



Tr DT" Tr Ei Tr E2 



N 



{Tr BB Tt C Ti C + 2Tv ADTr EiTv E2) + ... 



NTrATrA + 2TrATrBTrC + 



{TiBTiBTiCTTC + 2TrATiDTiEiTiE2) + ... . (4.19) 



For a meson multiplct other than the Goldstonc bosons, ml ^ O (1) and 
^i/m2 ~ (1/-/V). Using this in eq. (4.19) shows that the terms induced by 
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integrating out a heavy meson nonet have the same A^-counting as those 
in the original Lagrangian eq. (4.16), as one might expect. 

The singlet meson plays an important role in reproducing the correct N- 
counting, as there are non-trivial cancellations between singlet and octet 
meson exchange. It is inconsistent to include the octet mesons but not 
the singlet meson in the large N limit. Neglecting the singlet meson is 
equivalent to letting mg oo, so that Si/„i2 ~ — l/m| and is of order one. 
In this case, the ^i/^^ terms in eq. (4.19) violates the A'^-counting rules by 
one power of N. 

It is also inconsistent to use large counting rules, and treat the rj' 
as heavy. Integrating out the r]' is equivalent to retaining only the singlet 
meson contributions in eq. (4.19), which can be done by letting mg — > oo. 
In this case = 1/mg = l/m^, . The rj' exchange terms violate A'^- 

counting by two powers of A^, if one uses m^, ~ O {l/N). It is inconsistent 
to integrate out the rj' , and at the same time assume that m^, ~ C (1/A^), 
since a light particle is being integrated out of the effective Lagrangian. 
This is the origin of L7 ~ O (N^) ■ Retaining the rj' in the effective La- 
grangian gives Lt ^ O (1) [17]. 

5. Baryons 

Baryons are color singlet hadrons made up of quarks. The SU{N) invariant 
e-symbol has A^ indices, so a baryon is an A^-quark state, 

A baryon can be thought of as containing A^ quarks, one of each color, 
since all the indices on the e-symbol must be different for it to be non- 
zero. Quarks obey Fermi statistics, and the e-symbol is antisymmetric in 
color, so the baryon must be completely symmetric in the other quantum 
numbers such as spin and flavor. 

The number of quarks in a baryon grows with N, so one might think 
that large A^ baryons have little to do with baryons for A^ = 3. However, 
we will soon see that for baryons, as for mesons, the expansion parameter 
is 1/A'', and that one can compute baryonic properties in a systematic 
semiclassical expansion in l/N. The results are in good agreement with 
the experimental data, and provide information on the spin-flavor structure 
of baryons. We will be able to compute baryon properties such as masses, 
magnetic moments and axial couplings. The 1/7V expansion provides some 
deep connections between QCD and two popular models, the quark model, 
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t, 



(a) 



(b) 



Fig. 29. A baryon interaction and the corresponding connected component. 



and the Skyrme model, which provide a good phenomenological description 
of baryonic properties. 

5.1. N -Counting Rules for Baryons 

The A'^-counting rules for baryon graphs can be derived using our previous 
results for meson graphs. Draw the incoming baryon as A^-quarks with 
colors arranged in order, -N. The colors of the outgoing quark lines 
arc then a permutation of 1 • • • TV. It is convenient to derive the iV-counting 
rules for connected graphs. For this purpose, the incoming and outgoing 
quark lines are to be treated as ending on independent vertices, so that the 
connected piece of fig. 29(a) is fig. 29(b). A connected piece that contains 
n quark lines will be referred to as an n-body interaction. The colors on 
the outgoing quarks in an n-body interaction are a permutation of the 
colors on the incoming quarks, and the colors arc distinct. Each outgoing 
line can be identified with an incoming line of the same color in a unique 
way. One can now relate connected graphs for baryons interactions with 
planar diagrams with a single quark loop. The leading in N diagrams for 
the n-body interaction are given by taking a planar diagram with a single 
quark loop, cutting the loop in n places, and setting the color on each cut 
line to equal the color of one of the incoming (or outgoing) quarks. For 
example, the interaction in fig. 29(b) is given by cutting fig. 10 once at 
each of the three fermion lines. Planar meson diagrams contain a single 
closed quark loop as the outer edge of the diagram. Baryon n-body graphs 
obtained from cutting the quark loop require that one twist the quark lines 
to orient them with their arrows pointing in the same direction, and do 
not necessarily look planar when drawn on a sheet of paper. For example, 
fig. 30 is a "planar" diagram for a two-body interaction. Baryon graphs in 
the double-line notation can have color index lines crossing each other due 
to the fermion line twists. 

The relationship between meson and baryon graphs immediately gives 
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Fig. 30. An example of a "planar" two-body baryon graph. 



-1- 



Fig. 31. An example of a disconnected baryon interaction. 



US the A''-counting rules for an n-body interaction in baryons: an n-body 
interaction is of order N^~"', since planar quark diagrams are of order 
A'', and n index sums over quark colors have been eliminated by cutting 
n fcrmion lines. Baryons contain A'^ quarks, so n-body interactions are 
equally important for any n. n-body interactions are of order N^~'^, but 
there are 0{N") ways of choosing n-quarks from a A''-quark baryon. Thus 
the net effect of n-body interactions is of order N. 

Diagrams with two disconnected pieces, such as fig. 31, have a net effect 
of order 7V^, those with three disconnected pieces are of order N^, and 
so on. This is easy to understand. The baryon mass Mb is of order N, 
since it contains N quarks. The diagrams are an expansion of the baryon 
propagator, and sum to give 

Diagrams with a single connected component produce the order t term (the 
interaction can take place at any time) and are order TV ~ Mb, those with 
two connected components give the order term (each connected compo- 
nent can take place at any time) and are order N"^ ^ ^b^ Baryon 
interactions in the large limit are best studied in terms of connected dia- 
grams, and the diagrammatic methods are the same as used in many-body 
theory. 

Interactions of quarks in a baryon can be described by a non-relativistic 
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Fig. 32. Baryon matrix elements of a one-body operator such as qq, and a two-body 
operator such asqqqq. 



Hamiltonian if the quarks are very heavy. The Hamiltonian has the form 

H = Nm + y ^ + ^Yv(xi-Xj) 
^ 2m N ^ ^ ^' 

+ X] y {xi- Xj,Xi- Xk) + ■ . . . (5.1) 

Each term contributes 0{N) to the total energy. The interaction terms 
in the Hamiltonian eq. (5.1) are the sum of many small contributions, so 
fluctuations are small, and each quark can be considered to move in an 
average background potential. Consequently, the Hartree approximation is 
exact in the large N limit. The ground state wavefunction can be written 
as 

N 

iPo{xi,...,xn) = '[[(t>o {Xi), 

where the positions of the quarks. The spatial wavefunction 0o (x) 

is TV-indcpcndent, so the baryon size is fixed in the N ^ oo limit. The first 
excited state wavefunction is 

tpi{xi,...,XN) = -7=^(t>i{xk) Yl ^o{xi). (5.2) 

Further details about this approach can be found in refs. [3,5]. 

The iV-counting rules can be extended to baryon matrix elements of 
color singlet operators. Consider a one-body operator, such as qq. The 
baryon matrix clement {B\qq \B) has N terms, since the operator can be 
inserted on any of the quark lines (see fig. 32(a)). The baryon matrix 
element is therefore < 0{N). One obtains an inequality because there can 
be cancellations between the possible insertions. These cancellations will 
be crucial in unraveling the structure of baryons. Similarly, a two-body 
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Fig. 33. Diagrams for baryon-meson couplings. 




(a) (b) 
Fig. 34. Diagrams for baryon-meson scattering. 



operator such as qqqq has matrix element < 0{N'^), since tlicre are 0{N'^) 
ways of inserting the operator in a baryon (see fig. 32(b)). In general, an 
n-body operator has matrix elements < 0{N"). 

The baryon-meson coupling constant is < \/N. This can be seen from 
fig. 33, which shows the matrix element of a fcrmion bilinear in a baryon. 
There are N possible insertions of the fermion bilinear, so the matrix ele- 
ment is order N. The amplitude for a fermion bilinear to create a meson 
is the Z-factor, which is order \/iV, so the baryon-meson coupling con- 
stant is of order N/Z = Vn. The baryon-meson scattering amplitude is 
< 0{1). Two contributions to the scattering amplitude are shown in fig. 34. 
Figure 34(a) has N possible insertions of the fermion bilinear, and two me- 
son factors of each, so the net amplitude is < O (1)- The two 
bilinears must be inserted on the same quark line to conserve energy — the 
incoming meson injects energy into the quark line, which must be removed 
by the outgoing meson to give back the original baryon. If the bilinears 
are inserted on different quark lines, as in fig. 34(b), an additional gluon is 
needed to transfer energy between the two quark lines. The number of ways 
of choosing two quarks is 7V^, the meson factors are l/\fN each, and 
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Fig. 35. Diagrams for heavy meson + baryon — » heavy baryon. The heavy quark is 
represented by a double line. 



the two gluon couplings give each, so the total amplitude is again 

< 0(1). We have seen that the amplitudes baryon baryon + meson is of 
order ^/N , and baryon + meson baryon + meson is of order unity. One 
can similarly show that baryon + meson — > baryon + 2 mesons is of order 
1/Vn, etc. As in purely mesonic amplitudes, each additional meson gives 
a factor of 1/ a/ZV suppression. 

One can also look at the transition amplitudes for ground state baryon 
+ meson excited baryon, or equivalently, for B + Mq Bq, where Mq 
and Bq are mesons and baryons containing a single heavy quark. In both 
processes, one of the quarks in the final state is different from the others. In 
the transition amplitude diagram, the meson amplitude must be inserted 
on the quark line that is different, as in fig. 35; the meson operator either 
adds energy to the quark or converts it from a light quark to a heavy quark. 
Thus the combinatorial factor is unity, instead of A'^. The baryon with one 
excited (or heavy quark) has a wavefunction of the form eq. (5.2) in which 
one sums over the N possible quarks which can be different, and multiplies 
by a normalization factor of 1/\/N. This produces an additional factor of 
N X 1/VN, so the ^ M + B* amplitude is of order VN. This is 1/VN 
suppressed relative to the corresponding amplitude between ground state 
baryons. 

Baryon-baryon scattering amplitudes at fixed velocity are of order N. It 
is important to study baryon scattering at fixed velocity, rather than fixed 
momcntiim, because the baryon mass is of order TV. Working at fixed ve- 
locity avoids kinematic enhancements near threshold. The baryon-baryon 
scattering amplitude from diagrams such as fig. 36 has a combinatorial 
factor A^^ for the choice of two quarks, and (l/VlV)^ for the two gluon 
couplings, for an overall factor of A^. One could also consider a similar 
diagram without the exchanged gluon. Then the two quarks exchanged 
must have the same color, so the combinatorial factor is A^. The net result 
is that the baryon-baryon scattering amplitude at fixed velocity is of order 
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Fig. 36. Baryon-bciryon scattering. 



TV. Baryon-baryon scattering can be described by classical trajectories in 
the large N limit, since the particle masses are order N, and the scattering 
amplitude is also of order N . 

The processes considered so far all have an N dependence that is some 
power of N . There are also processes that are exponentially small in N , 
such as the cross-section for e^e~ — > BB. The amplitude to create a quark 
pair from the vacuum is some number a < 1. The baryon has N quarks, 
so the amplitude to create a BB pair is of order a^, and is exponentially 
suppressed in N. 

An important observation due to Wittcn is that all the A^-counting rules 
mentioned above arc the same as in a field theory with coupling constant 
I/a/Iv, where the mesons are fundamental fields and the baryon is a soliton. 

5.2. The Non-Relativistic Quark Model 

The non-relativistic quark model treats the baryon as made of three non- 

relativistic quarks boimd by a potential. The precise details of the potential 
will not be important for these lectures. All we need assume is that the 
ground state baryon is described by all three quarks in the same spatial 
wavefunction (t>{x). The wavefunction must then be completely symmetric 
in spin and flavor. In the case of three flavors, there are six possible quark 
states u'\,ul, d'\,dl, s'\,sl. The potential is assumed to be spin and 
flavor independent, so the non-relativistic quark model has an SU{6) spin- 
flavor symmetry under which these six states transform as the fundamental 
representation. The ground state baryons transform as the completely 
symmetric product of three 6's of SU{6), which is the 56 dimensional 
representation. 
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Fig. 37. Flavor SU{3) weight diagram for the decuplet baryons. The horizontal axis is 
I3, and the vertical axis is hypercharge. 



Three spin 1/2's added together can give spin 1/2 or spin 3/2, so the 
baryon 56 contains spin- 1/2 and spin-3/2 baryons. The spin-3/2 wavefunc- 
tion in the m = 3/2 state is simple, it is Iff, and is completely symmetric. 
The flavor wavefunction must also be completely symmetric. It can have 
the form uuu, {uud + udu + duu)/\/^, etc. The spin-3/2 baryons are the 
decuplet baryons, shown in fig. 37. 

The spin- 1/2 baryon wavefunctions are slightly more complicated. There 
are no spin- 1/2 baryons in which all three quarks are the same, for then 
the wavefunction would be completely symmetric in flavor, thus completely 
symmetric in spin, and so spin-3/2. Consider a spin-1/2 baryon in which 
two of the quarks are identical, such as uud. The spin wavefunction for 
the two identical quarks must be completely symmetric, so the total spin 
wavefunction of the baryon in an to = 1/2 state must have the form 

TTi+A(Ti + iT)T. 

The constant A = —1/2 can be determined by requiring that the raising 

operator J+ annihilates the state, since it is a j = 1/2, to = 1/2 state. 
Thus the wavefunction of the baryon can be written as 

uwd [2 TTi - TiT - in] . (5.3) 

Actually, one should add cyclic permutations (and divide by 1/ VS) to en- 
sure that the wavefunction is completely symmetric. However, for most 
calculations, one can work just as well with eq. (5.3). We have determined 
the wavefunctions of six of the octet baryons. The remaining two states 
have three different quarks, uds. The S° is constructed using the combina- 
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Fig. 38. Flavor SU(3) weight diagram for the baryon octet. The horizontal axis is I3, 
and the vertical axis is hypercharge. 

tion {uds + dus) / \f2. This symmetrizes the wavefunction in the first two 
flavors, and so one constructs the wavefunction as in eq. (5.3), 

^{uds-rdus) [2TTi-TiT-iTT]. 

The A state is constructed by antisymmetrizing uds in the first two fiavors, 
{uds-dus)l\pi. The spin- wavefunction must also be antisymmetric in the 
first two flavors, so the A wavefunction is 

\ {uds - dus) [tit - in] , 

which can be abbreviated to 

-^uds [Tit- ITT] • 

The entire spin of the A is carried by the s-quark. The spin-1/2 octet is 
shown in flg. 38. The spin-1/2 octet and spin-3/2 decuplet together make 
up the 56 of SU{6). The permutation symmetry properties of the baryons 
under spin SU{2) and flavor SU{3) are: 

8 = (□, 10 = (cm, cm) . 

It is straightforward to compute baryon properties in the non-relativistic 
quark model. The axial coupling constant qa of the proton is given by the 
matrix element (p| q^^^^T^ \p). In the non-relativistic limit, this reduces 
to the matrix element of q^a^T^q, an operator which is +1 on u t and d J,, 
— 1 on u J, and d t and zero on strange quarks. The proton matrix element 
is 

5^ = ^[l + l + l] + i[l-l-l] + l[-l + l-l] = | 

where the first term is (2/\/6)^ times the matrix elements for u'\,u'\, d \., 
etc. obtained using the wavefunction eq. (5.3). 
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Fig. 39. Flavor SU{3) weight diagram for the spin-1/2 baryons for N colors. The 
numbers represent the degeneracy of each weight. The long edge of the diagram has 
{N + l)/2 states. 



Similarly, one can compute the magnetic moments of the baryons. The 
magnetic moment operator is /Lkt^, where ^ is the magnetic moment of the 
quark, and is the spin operator. The quarks have magnetic moments 
A*M) A*d and jjg. If the u and d quarks are degenerate in mass, /U„ = —l^^. 
The computation of the baryon magnetic moments in the non-relativistic 
quark model is left as an exercise. 

One can repeat the entire non-relativistic quark model analysis for A'^ 
colors. It is convenient to choose A'' = 2m+l, so that N is always odd. The 
baryons form the completely symmetric tensor of SU (6) with N indices, 

iim---m (5.4) 

This decomposes under SU (2)spin x 5C/(3)flavor as a tower of representations 
with J = 1/2, J = 3/2, . . . J = N/2, 

(i.ra-p) (f.ffl-^)-(f.™-™) (5.5) 

The SU{3) representations are complicated for arbitrary values of N. For 

example the SU{3) weight diagram of the J = 1/2 baryons is shown in 
fi-g. 39. The flavor representations simplify for two light flavors, where the 
states are 7 = J = 1/2, ...I = J = N/2. 

Baryon transformation properties under spin and flavor are N dependent 
for Np > 3, unlike for mesons. To apply the 1/N expansion to baryons it 
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is convenient to make some identification between the arbitrary N baryons 
and the = 3 baryons. The proton state for arbitrary N will be taken 
to be the strangeness zero 7 = J = 1/2 state. It contains (m + 1) u 
quarks and m d quarks. The spin of the u quarks must be (m + l)/2, 
since the spin wavefunction has to be completely symmetric. Similarly, the 
spin Jd, of the d quarks must be m/2. The proton is then the J = 1/2 
state made by combining ,/„ and Jd to form spin-1/2. This is sufficient 
information to compute many of the proton properties as a function of 
N [31]. For example, one can compute qa, which is the matrix element of 
(j^T^ = 2{Jl- Jl). By the Wigner-Eckart theorem 

(p|Ju \p) = K (p| Jb), 

(j>\ Jd \p) = Xd (p\J\p), 
so that 

> J ■ Ju > J ■ Jd 

■'^u — — ^9 — ) -^d — 



j2 , j2 • 

Using J = J„ + Jrf, one finds that 

2J.J„=J2+J2 -J2, 2J.Jrf = j2+j2-j2, 

so that 

_ + 5 A^-l 

A« — T , Ad — • 

6 6 

In the large N proton, u-quarks tend to have spin up, and d-quarks tend 
to have spin down. The axial coupling is 

9a = Au — Ad = — ^ — , 

which reduces to 5/3 when = 3. gA ^ 0{N) in the large limit, 
which is consistent with the iV-counting rules for the matrix element of a 
one-quark operator. 

Problem 5.1 (Baryon Magnetic Moments) 

Show the following: 

(i) The magnetic moment of a spin-1/2 baryon qaQaQb with two identical 

quarks is 4^a/3 — l-ib/S 

(a) The magnetic moment of the A is fig 

(Hi) The magnetic moment of the E° is the average of the S+ and S~ mag- 
netic moments. 
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(iv) The magnetic moments of the spin-3/2 baryons is the sum of the mo- 
ments of the constituent quarks. 

(v) Find the A+ transition magnetic moment. 

6. Spin-Flavor Symmetry for Baryons 
6.1. Consistency Conditions 

The large A'^ counting rules for baryons imply some highly non-trivial con- 
straints among baryon couplings. The simplest to derive are relations 
among pion-baryon couplings, or equivalently, baryon axial current ma- 
trix elements. Related results also hold for p-baryon couplings, etc. and 
are discussed later. To derive the axial current relations, consider pion- 
nucleon scattering at fixed energy in the N ^ oc limit. The argument is 
simplest in the chiral limit where the pion is massless, but this assumption 
is not necessary. The two assumptions required are that the baryon mass 
and qa are both of order N. One expects the baryon mass to be propor- 
tional to N, since it contains N quarks. We have seen that the A^-counting 
rules imply that gA is order A^, unless there is a cancellation among the 
leading terms. In the non-relativistic quark model, gA = {N + 2)/3, so 
such a cancellation does not occur. It is reasonable that gA is of order N 
in QCD, even though it need not have the value {N + 2)/3. 
The pion-nucleon vertex is 

^(B|g7^5T»(z|B), 

and is of order VTV, since gA ^ N and /tt ^ \/N. Recoil effects are of order 
1/A'^, since the baryon mass is order N and the pion energy is order one, 
and can be neglected. This allows one to simplify the expression for the 
nucleon axial current. The time component of the axial current between 
two nucleons at rest vanishes. The space components of the axial current 
between nucleons at rest can be written as 

{B\ ^iSsT^V \B) = gN {B\ X» \B) , (6.1) 

where {B\ X^" \B) and g are of order one. The coupling constant g has been 

factored out so that the normalization of X*" can be chosen to simplify 
future expressions. X^°' is an operator (or 4x4 matrix) defined on nucleon 
states p '\, p i, n '\, n I, and X^" has a finite N ^ oo limit. 

The leading contribution to pion-nucleon scattering is from the pole 
graphs in fig. 40, which contribute at order E provided the intermediate 
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Fig. 40. Pion-nucleon scattering diagrams of order E. The third diagram is l/AT^ 
suppressed in the large N limit. 



state is degenerate with the initial and final states. Otherwise, the pole 
graph contribution is of order E'^. In the large N limit, the pole graphs are 
of order N , since each pion-nucleon vertex is of order ^/N. There is also a 
direct two-pion-nucleon coupling that contributes at order E, which is of 
order 1/N in the large N limit and can be neglected. 

The pion-nucleon scattering amplitude for Tr°'{q)+B{k) — > Tr^{q')+B{k') 
from the pole graphs is 



N' 



1, 



J TT 



1 j^j5^ia _ 1 Jfiaj^jfc 



(6.2) 



where the amplitude is written in matrix form, with the matrix labels de- 
noting the spin and flavor quantum numbers of the initial and final nucle- 
ons. Both initial and final miclcons are on-shell, so ~ q'^ . The product 
of the X matrices in eq. (6.2) sums over the possible spins and isospins 
of the intermediate nucleon. Since /^r ^ ^/N, the overall amplitude is of 
order A^, which violates unitarity at fixed energy, and also contradicts the 
large N counting rules of Witten. Thus large N QCD with aJ = J = 1/2 
nucleon multiplet interacting with a pion is an inconsistent field theory. 
There must be other states that cancel the order N amplitude in eq. (6.2) 
so that the total amplitude is order one, consistent with unitarity. One 
can then generalize X"* to be an operator on this degenerate set of baryon 
states, with matrix elements equal to the corresponding axial current ma- 
trix elements. With this generalization, the form of eq. (6.2) is unchanged, 
and we obtain the first consistency condition for baryons [32], 

=0 . (6.3) 

This consistency condition implies that the baryon axial currents are rep- 
resented by a set of operators which commute in the large N limit, 
a result also derived by Gervais and Sakita. [33]. There are additional 
commutation relations, 

[j\X^^]^ie,jkX''\ (6.4) 
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since has spin one and isospin one. 

The algebra in eqs. (6.3) and (6.4) is a contracted SU{2Np) algebra, 
where Np is the number of quark flavors. To see this, consider the algebra 
of operators in the non-relativistic quark model, which has an SU{2Nf) 
symmetry. The operators are 

where J' is the spin, T" is the flavor generator, and G'" are the spin-flavor 
generators. The commutation relations involving G"* are 

[G*", G^''] = -^j^ eijk6ab j'' + fabcSij ^'^ + 2 ^ijkdabc G''" , 

[T^G*^] =z/„fccG^^ (6.5) 
The algebra for large N baryons in QCD is given by taking the limit 

X'«= lim — . (6.6) 

The SU{2Nf) commutation relations eq. (6.5) turn into the commutation 
relations cqs. (6.3 6.4) in the large N limit. The limiting process eq. (6.6) 
is known as a Lie algebra contraction. 

We have just proved that the large N limit of QCD has a contracted 
SU {2Np) symmetry in the baryon sector. The unitary irreducible repre- 
sentations of the contracted Lie algebra can be obtained using the theory 
of induced representations, and can be shown to be infinite dimensional. 
The simplest irreducible representation for two flavors is a tower of states 
with I = J = 1/2, 3/2, etc. which is the set of states of the Skyrme model, 
or the large N non-relativistic quark model. The irreducible representa- 
tions for three flavors are more complicated. The 1/7V expansion allows one 
to compute baryonic quantities using SU(2Np) symmetry in the N ^ co 
limit. The 1/N corrections allow one to systematically study the form of 
SU{2Nf) symmetry breaking at finite A^. 

The pion-baryon coupling matrix X*" can be completely determined 
(up to an overall normalization g), since it is a generator of the SU{2Np)c 
algebra. It is easy to show that the large QCD predictions for the 
pion-baryon coupling ratios are the same as those obtained in the Skyrme 
model or non-relativistic quark model [34] in the N ^ oo limit, because 
both these models also have a contracted SU{2Np) symmetry in this limit. 
In the Skyrme model, the axial current in the N ^ oo limit is oc 
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Tr AT'^A^^T"- , where A is the Skyrmion coUective coordinate. The X's 
commute (and so form part of an SU [2Nf)c algebra), since A is a c-number. 
We have already seen how the quark model algebra reduces to SU{2Np)c 
in the large N limit. While we have shown that SU {2Np)c is a symmetry of 
QCD in the large N limit, we have not shown that SU{2Np) is a symmetry 
of QCD for finite A''. There is no reason to believe this is the case, so 
the SU{2Np) symmetry of the quark model is not a symmetry of QCD. 
Nevertheless, many results obtained in the quark model will be rederived 
in QCD using the SU{2Np)c symmetry that is exact when N ^ oo. 

It is useful to have an explicit representation of the N ^ oo SU{2Nf)c 
algebra to compute baryon properties in the expansion. There are 
two natural choices discussed above, the quark model (acting on N quark 
baryons) , or the Skyrme model. The two give equivalent results for physical 
quantities. An operator O (such as X*") in the quark representation can 
be written as a 1/A'' expansion of Skyrme operators, and vice- versa, 

Oquark = Og^^^ + "^Oskyrme + i^''^) 

A typical expression for a physical quantity is 
aoOo + ^0, + ^0, + ... 

where are coefhcients, and Oj are operators (such as or J*, etc.) 
in either the quark or Skyrme representation. The coefficients a, are not 
the same in the two representations, and the relation eq. (6.7) can be 
used to relate them. Clearly, any difference between the representations 
is of higher order in 1 /N than the terms retained in the calculation. This 
is similar to the scheme dependence of physical quantities in Feynman 
diagram perturbation theory. 

The 1/N expansion provides some deep connections between the quark 
model, the Skyrme model, and QCD. I do not have time in these lectures 
to discuss both the quark and Skyrme models, and the connection between 
them. Some of these results can be found in refs. [34-36]. In the rest of 
these lectures, I will compare the large N results mostly with the quark 
model. 



Problem 6.1 {SU{2Np) Commutation Relations) 

Derive the SU{2Np) commutation relations eq. (6.5) using the quark op- 
erator relation [q,q^] =1. 
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Fig. 41. The leading contribution to nN — > irnN in the large N and chiral limits. 



6.2. 1/N Corrections 

What makes the l/N expansion for baryons interesting is that it is possible 
to compute corrections. This allows one to compute results for the 
physically relevant case N = 3, rather than for the strict N = oo limit, 
which is only of formal interest. 

The corrections to the axial couplings X*" are determined by con- 
sidering the scattering process tt" + B — > tt** + tt'^ + B at low energies. The 
baryon pole graphs that contribute in the large N limit are shown in fig. 41. 
The axial coupling X^°- can be expanded in a series in 1/A'', 



(6.8) 



where Xq" satisfies eq. (6.3). An explicit calculation shows that the ampli- 
tude for pion-nucleon scattering from the diagrams in fig. 40 is proportional 
to 

^3/2 [^^a^ [X1\X'"']] 

times kinematic factors, and violates unitarity unless the double commu- 
tator vanishes at least as fast as N~^/^, so that the amplitude is at most 
of order one. In fact, one expects that the double commutator is of order 
from the large TV counting rule that the amplitude is order 1/\/N. 
Substituting cq. (6.8) into the constraint gives the consistency condition 



' 



X 



\rkc 

^0 



X, 



' 



T^jfc ykc 



0, 



(6.9) 



usmg 



via yjt 



from eq. (6.3). For two flavors, the only solution to 

eq. (6.9) is that XJ" is proportional to Xg" [32]. This can be verified by 
an explicit computation of Xl°- as shown in section 8 [32,35]. In the rest 
of this section, I will state the solutions of the various large N consistency 
conditions. The solutions can be obtained using the methods discussed in 
detail in section 6.3. 
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At order 1/N, the baryons in an irreducible representation of the con- 
tracted SU{2Np) Lie algebra are no longer degenerate, but are split by 
an order 1/N mass term AM. The intermediate baryon propagator in 
eq. (6.2) should be replaced by 1/{E — AM). The energy E of the pion is 
order one, whereas AM is of order so the propagator can be expanded 
to order 1/N as 

1 1 AM 



Including the 1/N corrections to the propagator does not affect the deriva- 
tion of eq. (6.3), as the two terms in eq. (6.10) have different energy de- 
pendences. The first term leads to the consistency condition eq. (6.3) and 
the second gives the consistency condition on the baryon masses [37,35] , 

[X'", AM]]] = 0. (6.11) 

This constraint can be used to obtain the 1/A'' corrections to the baryon 
masses. The constraint eq. (6.11) is equivalent to a simpler constraint 
obtained by Jenkins using chiral perturbation theory [37] 

[X'"", [X"',AM]] = constant. (6.12) 

The solution of eq. (6.11) or (6.12) is that the baryon mass splitting AM 
must be proportional to J'^/N = j{j + 1)/N, where j is the spin of the 
baryon [37]. This is precisely the form of the baryon mass splitting in the 
Skyrme model [38] or non-relativistic quark model. 

The structure of the 1/A'' corrections shows that the expansion param- 
eter is J/N, where ,/ is the spin of the baryon. For example, the baryon 
mass spectrum including the J^/N mass splitting has the form shown in 
fig. 42. The correction terms are only small near the bottom of the (infinite) 
baryon tower. For this reason, the 1/A^ expansion will only be considered 
for baryons with spin J held fixed as ^ oo. 



6.3. Solution of Consistency Conditions 

The solution of the large N consistency conditions is the key to under- 
standing the structure of baryons. The answer can be given quite sim- 
ply [35,36,39,40]. A quantity Q that is of order N^ in the large A^ limit can 
be expanded in terms of the basic operators of the quark representation, 
G'", J* and T" as 
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Fig. 42. The baryon mass spectrum including the J'^/N term. The top of the tower is 
j = N/2, and the bottom is j = 1/2. 



where T' is a polynomial in its arguments, with coefficients that have an 
expansion in For example, we will soon see that the baryon mass in 

the flavor limit has the expansion 



Mb 
N 



(6.14) 



where the ai{l/N) are unknown expansion coefficients. It is important 
to remember that no assumption is being made about the validity of the 
non-relativistic quark model. The operators G"°, and are used as 
a way to do the SU{2Nf) group-theoretic computations. The expansion 
eq. (6.13) is true irrespective of the quark mass. For very heavy quarks, the 
non-rclativistic Hartrec picture is valid, and one can see that eq. (6.13) is 
consistent with the diagrammatic analysis oi l/N factors due to gluon ex- 
change [39,40]. Two-body operators on the right-hand side of eq. (6.13) are 
generated by an insertion of a one-body operator plus one-gluon exchange, 
three-body operators by a one-body operator plus two gluon exchanges, 
and so on, as shown in fig. 43. Each gluon exchange brings with it a fac- 
tor of l/N from the two couplings, which reproduces the A'^-counting in 
eq. (6.13). What is non-trivial is that the A^-counting rules also hold if the 
quarks are massless. 

At first sight, all terms in the expansion eq. (6.13) arc equally important, 
since an r-body operator in the expansion of V has a coefficient of order 
1/N^, and matrix elements of r-body operators are of order A''''. Including 
all terms in eq. (6.13) is cqtuvalcnt to saying that all possible SU{2Nf) 
representations are equally important, so that there is no predictive power. 
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Fig. 43. Two-body and three-body operator contributions to the baryon matrix elements 
of a QCD operator. 



Table 2 

51/(6) Identities. Some of the identities need to be projected onto a given spin 
or flavor channel, which is given in parentheses. 



2{J\J'}+3{T'',T''} + 12{G''',G»"} = 5N{N + 6) 



^abc^QiaQibj ^ 2 yi^Qicj + 1 d^i-c {T", T''} = | (AT -f 3) T= 



-12 {G'^.G'"} +27 {r",T"}-32 {J\J'}=0 




^abc {Gia^QiiJ + 1 ^a6c j^J^a J'b^^ _ 1^ j^ji^Qic^^Q 




4 {G''',G'''} = {r°,r''} 


(27) 


gijfc G^^J- JcLbc ^J^^^ G^^} 




2 fj^bc 1 'J^a G^^l (1^3^ Jabc | ^ia (^jb^ 




^ijk ^Q''°-^Qjb'^ _ jacg^bch !^X^,G'''^} 


(TO-HIO) 


3 {G'^.G^"} = {J% JJ } 


(J = 2) 


3 <i°''= {G*",G^*'} = {J%G^=} 


(J = 2) 



There are operator identities which allow one to simplify the general ex- 
pansion cq. (6.13). and drop certain terms as sublcading in 1/N. It is these 
identities which allow one to use the 1/N expansion to make non-trivial 
predictions for baryons in QCD. The complete set of operator identities 
was derived in rcf. [36]. They are listed in Tables 2 and 3 for the case of 
three and two flavors, respectively. 
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Tabic 3 

5i7(4) Identities. Some of the identities need to be projected onto a 
given spin or flavor channel, which is given in the second column. The 
last column gives the transformation properties of the identities under 
SU{2) X 5C/(2). 



{J*, J*} + {7<',7°} +4 {G*",©'"} = |Ar(iV + 4) 




(0,0) 


2{J',G"'} ={N + 2) 




(0,1) 


2{/",G'"} = {N + 2) r 


(1,0) 




i {j'',!"} -€^^^e'^^''{G^'^,Gi^} = (N + 2) G^" 




(1,1) 


{7°,7»} - {J\J^] = 


(0,0) 




4 {G'^.G'*"} = {/°,/^} 


(7 = 2) 


(0,2) 


fjjfc 1 J'jG^'^} = e"*"^ {/'',G''''} 




(1,1) 


4 {G*",G-'"} = {J% J-J } 


(J = 2) 


(2,0) 



The operator identities allow one to eliminate certain operator combi- 
nations in the general expansion eq. (6.13). For example, for three flavors, 

{T",G"} = \{N + 2,)r 

can be used to rewrite the two-body operator on the left-hand side in terms 
of the one-body operator on the right-hand side. The two-body operator 
matrix element is of order N'^^ and the one-body operator matrix element 
is of order N ^ which is consistent with the coefficient of proportionality 
being of order N . The operator identities have been written in terms of 
anticommutators, since they are hermitian, and since commutators can be 
simplified using the SU{2Np) commutation relations eq. (6.5). A study of 
the identities leads to the following reduction rules: 

Operator Reduction Rule (three flavors): All operator products in 
which two flavor indices are contracted using S"^^, d"-'"^ or /"''^, or two spin 
indices on G's are contracted using 6"^^ or e*-'*^ can be eliminated. 
Operator Reduction Rule (two flavors): All operators in which two 
spin or isospin indices are contracted with a ^ or e-symbol can be eliminated, 
with the exception of J^. 
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The inclusion of J^, but not in the set of independent operators does 
not break the symmetry between spin and isospin, because of the identity 

j2 ^ j2 * 

As an apphcation of the operator reduction rule, consider the baryon 
masses in the SU{3) symmetry limit. The general form of the baryon mass 
is given by eq. (6.13), 

where the terms in V must be spin- zero flavor singlets, since those arc the 
quantum numbers of M. Thus all terms in V are obtained by contracting 
the spin and flavor indices on G™, J' and T" using spin and flavor invariant 
tensors, such as 5°^'' or The operator reduction rule tells us that all 

terms involving or G*" can be eliminated, so that one is left with an 
expansion in J*. Rotational invariance implies that the expansion is in J^, 
and so has the form eq. (6.14), and leads to the spectrum in fig. 42. 

Equation (6.14) shows that the baryon expansion parameter is J/N. One 
can compute baryon properties in a systematic expansion in 1/A'' provided 
one looks at states with J fixed as — !■ oo. Generically, J is a one-body 
operator and can have matrix elements of order A^. The 1/N expansion 
is only valid for states in which J is of order unity, i.e. for states at the 
bottom of the baryon tower in fig. 42. In these states, there is a cancellation 
between the A^ possible insertions of the one-body operator J on the N 
quark lines. Generically, matrix elements of G"* and T" are 0{N), so 
that G'^°'/N and T"-/N arc 0(1). In the case of two flavors, there is an 
additional simphfication that /" is C (1), so that 1°- /N is 1/A^ suppressed. 

Problem 6.2 

Prove the identity 

AJ'G'" = {N + 2)1°' 
for two quark flavors. 

Solution: It is convenient to write G*" = (1/4) ''■q'''q; where the sum 
is over the different quark lines in the baryon, each of which has a different 
color. Then 

E^i ^o, I \ ^ „i „i 

* 

For two flavors, the isospin T" will also be denoted by I"^. 



66 



A. V. Manohar 



The term with a = (3 can be written as 

a=/3 a 

= 6r. 

The terms with a ^ (3 can he simphhed using the identity 

where Sexch{ce,P) is the spin-exchange operator that exchanges the spins 
of the two quarks, to give 

E ^ E [2^exch(a, /3) - 1] 

The final state baryon is completely symmetric in spin (K" flavor, so one 
can replace Sexch{ce,P) by the flavor exchange operator Fexch(Q;, /?). The 
identity 

[rS [rlYd = - m ^ 2Fe.ch(a,/3) - 1 
allows one to rewrite this as 

Then 

implies that 

E^w=^a('^r+^^s«'>^^') 

= E ^« + E ^^aahT'cT-"^ - E i^gahT^T^ 
a a j3 

a 

= 2{N - 1)P + Aiegahl'l'' + 2egaheghrr 
= 2{N - 1)P + 2iegah [/^ I''] + 2egaheghrr 
= 2{N-1)P 
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Combining the pieces gives the desired identity. More examples of this kind 
can be found in refs. [36,39,40]. 



7. Masses with SU(3) Breaking 

The baryon masses can be computed in a systematic expansion in powers 
of 1/A'' and SU{3) breaking. I will use a simplified version of the analysis 
of Jenkins and Lebed [41], by neglecting isospin breaking. This example is 
discussed in detail, to show how the formalism we have developed can be 
applied. 

There are eight isospin-avcragod baryon masses for the N, A, S, S, A, 
S*, S* and i7. The general form of the SU{3) singlet mass term has already 
been worked out in eq. (6.13). At first order in SU{3) breaking, the baryon 
mass term transforms as an SU{3) octet. The most general spin-zero SU{3) 
octet is a polynomial in J*, and G'" with one free flavor index set to 

8. All operators with contracted flavor indices can be eliminated using the 
operator reduction rule, so one is left with terms with a single or G*^, 
and powers of J', 

+e6,^+e64^^ + ... (7.1) 

where e is a measure of SU (3) breaking. SU (3) breaking is due to the 
strange quark mass mg. Chiral perturbation theory shows that there are 
non-analytic contributions to the baryon mass, such as TOs^^ and m!^ log m^. 
The structure of SU (3) breaking as a function of m., is highly non-trivial, 
but goes beyond the scope of these introductory lectures. For the purposes 
of the analysis here, I will treat e as a small parameter of order SU{3) 
breaking. 

At second order in SU (3) breaking, one can get a tensor with two free 
flavor indices set to 8, and at third order, with three free flavor indices set 
to 8. One can work out the general form of these terms, as in eq. (7.1). 
Eventually, the matrix elements will be taken between baryons for iV = 3. 
For this reason, one can stop the expansion of the mass term at three-body 
operators, since a baryon for A'^ = 3 contains three quarks. Operators with 
more than three quarks can be written as a linear combination of operators 
with three or fewer quarks. One simple way to see this is to write operators 
in normal ordered form, with all q^^s to the left of the q^s. Then r-body 
operators with r > 3 vanish on baryons containing three quarks. 
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The order terms up to three-body operators are 

e^C2^+e^C3^^ (7.2) 

and the term is 

^''^3^^. (7.3) 

Combining cqs. (6.13)-(7.3) and stopping at three-body operators gives the 
baryon mass formula 

j2 J*G'^ J^T^ J^SrpS 

M = ao + ai— + ebiT^ + e''2— ^ + e^3-^ + ^^^2-^ 

+^ + (7.4) 

which gives the eight baryon masses in terms of eight parameters. This 
must be the case, since the baryon masses are independent if one works 
to arbitrary order in and e. The non-trivial information contained in 
eq. (7.4) is the e and N dependence of the various terms. Equation (7.4) 
will be used to derive a hierarchy of baryon mass relations to a given order 
in e and To obtain these relations, one needs the matrix element of 

eq. (7.4) between baryon states. The matrix clement of 

1 \ / 1 




12 T« 



IS 

T' = ^iN-3Ns), (7.5) 



where Ng is the number of strange quarks. The operator G*^ is 

G** = (.r - 3 J!) , 



where J* is the quark spin, and 

^-''?(° » J' 

is the strange quark spin. This gives 

jiQiS ^ 1 / j2 _ 3 J . j \ 

v/12 ^ ^ 
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Table 4 

Matrix elements of the mass operators, eq. (7.4). The first section of the table 
lists the matrix elements of the basic operators N, Ng, J^, J| and 7^ . These 
axe used to compute the matrix elements of the remaining operators using 
eqs. (7.5)-(7.6). 





N 


A 


E 




A 


S* 






N 


3 


3 


3 


3 


3 


3 


3 


3 


Ns 





1 


1 


2 





1 


2 


3 


J2 


3/4 


3/4 


3/4 


3/4 


15/4 


15/4 


15/4 


15/4 


■/I 





3/4 


3/4 


2 





3/4 


2 


15/4 


/2 


3/4 





2 


3/4 


15/4 


2 


3/4 





1 


1 


1 


1 


1 


1 


1 


1 


1 




3/4 


3/4 


3/4 


3/4 


15/4 


15/4 


15/4 


15/4 




3 








-3 


3 





-3 


-6 


4^3 J' G'* 


3/2 


-3 


3 


-9/2 


15/2 





-15/2 


-15 




9/4 








-9/4 


45/4 





-45/4 


-45/2 


(2V3T^y 


9 








9 


9 





9 


36 


24T8j*G** 


9/2 








27/2 


45/2 





45/2 


90 




27 








-27 


27 





-27 


-216 



The total spin of the baryon is J = Jud + Js ■, where Jud is the spin of the 
u- and d-quarks, and Js is the spin of the strange quarks. One can write 
J"^ - 3J ■ Js = 3 J^^ - - 3Jg, and then use the identity J^^ = P, where 
/ is the isospin, to obtain 

jiQiS ^ 1 r^j2 _j2_ 3j2N .7 



All the baryons are eigenstates of A^^, P, and , so the matrix element 
of cq. (7.4) can be computed simply using eqs. (7.5) and (7.6). The matrix 
elements are listed in Table 4. 

Combining eq. (7.4) and Table 4 gives the baryon masses 

M = HA, (7.7) 

where 



M = thm niA tuy. ms ttia ttis* toh* ma ) , 
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Table 5 

The order in l/N and e of the eight 
mass operators. 





N 1 1/N 


i/Ar2 


1 


1 J2 




e 












e3 







is a row vector of baryon masses, 



A = 



f 


1 


1 


1 


1 


1 


1 


1 


1 


\ 




3/4 


3/4 


3/4 


3/4 


15/4 


15/4 


15/4 


15/4 






3 








-3 


3 





-3 


-6 






3/2 


-3 


3 


-9/2 


15/2 





-15/2 


-15 






9/4 








-9/4 


45/4 





-45/4 


-45/2 






9 








9 


9 





9 


36 






9/2 








27/2 


45/2 





45/2 


90 






27 








-27 


27 





-27 


-216 


/ 



is the array of matrix elements from Table 4, and 
H = [ Hi H2 H3 Hi H5 He Hr Hg ) 

~ V ^ 2\/3 4\/3iY 2V3N^ 12JV liN^ 24^3^ J 

is the row vector of coefficients. 

The classification of the mass operators in powers oil/N and e is shown 
in Tabic 5. There arc no relations if all the operators are retained in 
eq. (7.7). The most accmatc relation is obtained if one omits the operator 
T^T^T^ which contributes at order e'/N'^. This gives a baryon mass re- 
lation that has an error of order e^/N'^. The mass relation is obtained by 
writing eq. (7.7) as 

H = MA-'^. (7.8) 

Omitting T^T^T^ means that Hg, = 0, which gives the relation i?8 (on 
multiplying by 162 to eliminate fractional coefficients) 



A - 3S* + 3S* - n = O {e^/N^) , 



(7.9) 
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Table 6 

Baryon mass relations with the order in e and 1/N. The entries in the table are the 
coefficients of tum, etc. in the eight relations. The fractional error computed using 
the experimental values for the masses is listed in the last column. 





N 


A 


S 




A 


E* 




n 


Order 


Frac. Error 


Rl 


2 


1 


3 


2 


8 


6 


4 


2 


N 




R2 


-10 


-5 


-15 


-10 


16 


12 


8 


4 


1/N 


0.18 ±0.0004 


R3 


1 








-1 


4 





-2 


-2 


e 


0.27 ±0.0007 


R4: 


-7 


-2 


6 


3 


4 





-2 


-2 


e/N 


0.052 ±0.0003 


R5 


-2 


3 


-9 


8 


4 





-2 


-2 


e/Ar2 


0.011 ±0.0003 


R6 


2 


-3 


-1 


2 


16 


-20 


-8 


12 


e^/N 


0.0048 ±0.0004 


R7 


-14 


21 


7 


-14 


8 


-10 


-4 


6 




0.0017 ±0.0002 


R8 














1 


-3 


3 


-1 




0.0009 ±0.0003 



since ffs is O [e^/N^). The next most accurate relation is obtained by 
dropping T^T^T^ and T^J^G^^. This is equivalent to Hg = and Hj = 
in eq. (7.8). The relation Hs is eq. (7.9). The relation Hr = is 

{R7) : UN - 21A - 7S + 14S - 8A + lOE* - 4S* - 6f2 = O (eV-^^) • 

At the next step, one can drop cither the T^T^ operator, which gives a 
new relation with an error O (e^/N), or drop the J^T^ operator, which 
gives a relation with an error O {e/N^) ■ Dropping both gives no additional 
independent relation. One can then drop J'G*^ to get a O {e/N) relation, 
and then drop either or to get O (e) or O (l/N) relations. 

This procedure gives a hierarchy of mass relations Hi = 0. There are 
no free parameters in any of the relations, since all coefficients were fixed 
by group theory. It is convenient to write down relations that have defi- 
nite SU{6) and SU{2) x SU{3) transformation properties. These relations 
are orthogonal to each other, since different irreducible representations of 
a symmetry group are orthogonal. One can make the relations Hi = 
orthogonal using the Gram-Schmidt procedure starting with Hs = and 
working down to Hi =0, with respect to the metric (see problem 7.1) 

5 = diag( 1111^^11). (7.10) 

The entries in this matrix are the reciprocals of the number of baryon 
states, e.g. there are four nucleon states, p 1 , p l, n 1 , n two A states 
A t, A I, etc. The resulting relations are tabulated in Table 6. 

The mass relations can now be compared with experiment. The relations 
derived are homogeneous relations, and have no standard normalization, 
i.e. they can be multiplied by an arbitrary overall coefficient. It is therefore 
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best to compute the fractional error on each relation. Relations which hold 
to higher order in or e should work better, if and e are good 
expansion parameters. The fractional error is defined using the following 
procedure. Each relation is written in the form L — R, where both sides 
are linear combinations of masses with positive coefficients. For example, 
the O (e^/N'^) relation R8 is written as 

L = A + 3S* = 3S* + f2 = i?. 

One then computes the fractional accuracy of 
the relation, \L — R\ / {{L + R) /2) using the experimental values for the 
isospin averaged baryon masses (in MeV), 

rriN rriA ms mn 

938.91897 ±0.0002 1115.684 ± 0.006 1193.12 ±0.04 1318.11 ±0.31 

niA rriT,* nis* ma 

1231.3 ±1.1 1384.6 ±0.4 1533.4 ±0.3 1672.45 ±0.29 

Since the baryon mass is of order N, the fractional accuracy of -R8 is 
0(^e^/N^), since the denominator is order N. The fractional accuracy 
of the baryon mass relations is listed in Table 6 and plotted in fig. 44. The 
error bars on the points are from the experimental errors on the measured 
baryon masses. The points in fig. 44 have been plotted so that relations of 
the same order in e have the same symbol. The standard SU{3) analysis 
of baryon masses which might be familiar to some of you is equivalent to 
using the relations we have derived, but ignoring the powers of N. The 
Gell-Mann Okubo formula for the baryon octet, and the equal spacing rule 
for the baryon decuplet are linear combinations of the two relations and 
one e"^ relation we have obtained. Clearly, SU{'d) breaking alone is not 
the whole story, and including the factors of 1/iV provides a much better 
understanding of the data than SU{3) alone. It is obvious from the figure 
that the 1/N and e expansions explain the observed data. There is also 
clear evidence for the validity of the 1/N expansion. For example, the 
three order e relations are of the same order in SU(3) breaking, but differ- 
ent orders in 1/N, and the 1/7V suppression is obvious in the experimental 
data. One also gets new relations (such as Rl) which cannot be derived 
using SU{3), and work just as well as the relations derived using SU{3) 
symmetry. 

Problem 7.1 (Mass Relations) 

(i) Show that relations which transform as different SU{6) and SU{2) x 
SU (3) representations must be orthogonal with respect to the metric g 



Large N QCD 



73 



1 \ \ I r 



▲ 



- e 



o 



0.1 =- 



o 



0.01 



▲ 



o 



m 



0.001 ^ 



I 



2 



J \ \ I \ I I \ I 

4 6 
Mass Relation 




Fig. 44. Accuracy of the large N baryon mass relations from Table 6. The error bars 
are those on the experimentally measured baryon masses. The shape of the points gives 
their order in SU{3) breaking: O is C (1), A is O (e), □ is O (e^), and X is O (e^). The 
differences within the O (e) and O (e^) relations arc explained by including the factors 
of 1/7V. The accuracies of the baryon mass relations are l/Af^, e/N, e/N^, e/N^, e^/W-^, 
e^/N^ and e^/N^, respectively. 



in eq. (7.10). Hint: linear combinations of baryon states that transform 
as given weights of SU{6) representations are orthonormal combinations 
of the 56, i.e. when p '\ , p i, n '\ , n i, are used as the basis states. 

(a) Derive the mass relations in Table 6. 
8. Other Results for Baryons 

The procedure of the previous section can be used to analyze other baryon 
properties, such as magnetic moments and axial couplings [42-44]. The 
expansion has been applied to the nuclcon-nuclcon potential, and ex- 
plains the origin of Wigner supermultiplet symmetry in light nuclei [45,46]. 
Recently, there has been some interesting work on excited baryons using 
the expansion [47-49]. I do not have time to go over all these results 
in these lectures, and the reader is referred to the literature for details. In 
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this section, I will briefly discuss a few more results obtained using the 1/N 
expansion for baryons. 

The 1 /N analysis for two light flavors is much simpler than for three light 
flavors, since the baryon representations form representations of SU{4:)c 
symmetry, rather than SU{6)c- In the strangeness-zero sector, the baryon 

quantum numbers are / = J = 1/2, 3/2, These states will be referred 

to as a baryon tower. The S'C/(4) analysis can also be used in the case 
of more than two flavors. One can apply SU{4) spin-flavor symmetry to 
baryons in a given strangeness sector, so that the p and A are related to 
each other, the A, E, and E* are related to each other, and so on. What is 
more interesting is that one can relate different strangeness sectors to each 
other without assuming approximate SU{3) symmetry [37]. For example, 
requiring that the K + N ^ n + T, amplitude is unitary relates the pion 
couplings of the N and A to those of the A, E, and E*. It also constrains the 
form of the KNY, coupling [37,35]. Examples of relations obtained using 
large iV, but not assuming SU{3) symmetry are baryon mass relations such 
as [37,35,50] 

E*-E = S*-S-f 0(l/Ar2-)^ 

^^_^ = _l(0_S*-S*+A) + O(l/iV^), 

(S* - A) (fi - S*) = 2 (S* - S*) -h O (1/A^^) , 

i (E + 3E*) - A = ^ (A - TV) + O (l/iV^) , 

S^-S<3 = Sq-S'q + 0(1/A^'), 

1 (Eg + 3E^) - A= ^(A-iV) + {1/N^) , (8.1) 

and coupling constant relations such as 

gA = 9c,b + 0{l/N). (8.2) 

Equation (8.1) relates mass splittings of baryons of different strangeness, 

and also relates mass splittings of heavy quark baryons to mass splittings 
of baryons containing only light quarks. A more detailed analysis allows 
one to predict heavy baryon mass differences to high accuracy using the 
1/N expansion, by relating them to the known mass differences of the octet 
and decuplet baryons [50]. Equation (8.2) relates pion couplings of c and b 
baryons to the pion coupling of the nucleon, qa- This relation was originally 
derived using the Skyrme model in ref. [51]. One also obtains information 
on the Isgur-Wise function for heavy baryons using large N QCD [52]. 
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Table 7 

Values for the pion couplings from ref. [38]. The values of g^rNA and g^rNN 
axe obtained using the Skyrme model. The ratio g-nNii./9-nNN obtained 
using the Skyrme model is the same as that predicted using large N QCD 
to order l/AT^. 





Method 


Theory 


Experiment 


g-rvNA 


Skyrme Model 


13.2 


20.3 


g-rvNN 


Skyrme Model 


8.9 


13.5 


g-KNA/g-KNN 


Large N QCD/Skyrme Model 


1.5 


1.48 



In a given strangeness sector, one can predict the ratios of pion couplings 
to order 1 /N"^ . This result is simple to derive. Consider the operator expan- 
sion eq. (6.13) where the flavor group is now SU (2). The pion couplings are 
the same as the axial current matrix elements, by the Goldberger-Treiman 
relation. As in eq. (6.1), the only non-zero matrix element is that of the 
space component of the axial current. Thus the operator expansion is for 
spin one and isospin one. The expansion of the axial current has the form 

A'" = gG'" + + ..., (8.3) 

using the operator reduction rule for two flavors. The \/N expansion is 
applied to states with J of order unity. In the three-flavor case, matrix 
elements of G'" and T° can be of order N. However, in the two-flavor case, 
matrix elements of the isospin are of order unity, since = . The J*/" 
term is therefore a l/N"^ correction, so one can predict the ratios of axial 
current matrix elements (such as g^KNN / Q^KN^i.) with an accuracy of l/N"^. 
The Skyrme model is one particular representation of the SU{4:)c spin- 
flavor symmetry, so we have shown that the QCD predictions for the ratios 
of pion couplings, such as q-knn / Q-kNA or g^NN / QtvAA are equal to the pre- 
diction in the Skyrme model up to corrections of order 1/N'^. The Skyrme 
model results for the couplings are listed in Table 7. Clearly, the values 
of the individual couplings do not agree that well with the experimental 
data. Nevertheless, the ratio is in excellent agreement. Only the predic- 
tion for the ratio can be derived directly from QCD to O (l/iV^), and does 
not make any assumption about the validity of the Skyrme model. The 
overall scale of the couplings depends on the details of the Skyrme model 
Lagrangian, and is not a prediction of large N QCD. Similarly, other re- 
sults in the Skyrme model literature that are model independent (using the 
terminology of ref. [53]) can also be derived using large TV QCD. 

The axial coupling constant g in eq. (8.3) is the same for any baryon 
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Table 8 



Axial couplings extracted from 10 — i 


• 8 + TT 


decays. The last 


column gives the 


baryon 


strangeness. 






Decay 


9 


S 


A ^ AfTT 


1.8 





S* Ett 


1.5 


-1 


S* ^ Att 


1.5 


-1 




1.3 


-2 



tower, but it can be different for different towers, i.e. it can depend on 
tlie strangeness. One can prove that g is a, constant at leading order, and 
is linear in strangeness at order The value of g extracted from the 

different decuplet decays are given in Table 8, and clearly satisfy this result. 
One can also show that the F/D ratio for the baryon axial currents and 
magnetic moments is 2/3, with an error of order l/N"^. This is consistent 
with the experimental values of 0.58 and 0.59, respectively 

The I = J rule of Mattis and collaborators [54-56] can be derived in the 
expansion. This rule was originally derived using the Skyrme model, 
and is in good agreement with the experimental data. The I = J rule 
states that meson-baryon couplings satisfy I = J, where / is the isospin 
of the meson, and ,7 is the spin transfer at the vertex. For example, the 
p meson has 7=1, and its coupling to baryons is dominantly spin one, 
i.e. of magnetic moment type, or proportional to the F2 form factor. The 
CO has / = 0, and its coupling to baryons is dominantly spin zero, i.e. of 
charge type, or proportional to the Fi form factor. One can prove a slightly 
stronger form of the result using the 1/A'' expansion [36], the meson-baryon 

vertex is of order ( V ^ 1 

Finally, let me show two plots in which large N predictions have been 
compared with the experimental data. The first plot is fig. 45 for the 
baryon magnetic moment relations. The non-relativisitic quark model is 
known to provide a good description of the baryon magnetic moments at 
the 20% level. All the magnetic moment relations of the non-relativistic 
quark model can be derived using large N QCD to some order in 1/A''. It is 
clear from fig. 45 that the 1/iV expansion also explains why some relations 
work better than others; something that cannot be understood solely on 
the basis of the quark model. In addition, one gets one new relation (the 
filled circle) that that works just as well as the relations that are also valid 
in the quark model. The only relation that does not work as well as one 
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Fig. 45. Comparison of the baryon magnetic moment relations of ref. [42] with the 
experimental data. The circles are O {l/N) relations, the squares are O (l/AT^) relations, 
and the triangles are 0(e/N) relations. All the relations except the filled circle can 
also be derived in the non-relativistie quark model. The error bars are due to the 
experimental errors on the measured baryon magnetic moments. 



might expect is relation 7, which is an O (l/N"^) relation but is violated at 
the 25% level. This is a prediction for the A+ — » transition amplitude, 
and is a known problem for the quark model. 

The second plot is a comparison of the large N predictions for the 
nucleon-nucleon potential with the experimental data [46]. It is difficult 
to compare the predictions directly with nucleon scattering data. What is 
actually shown in fig. 46 is a comparison of the large TV predictions with 
coupling constants in the Nijmegen potential [57,58], which provides a good 
description of the experimental nucleon-nucleon scattering data. The 1/N 
expansion provides a satisfactory explanation of the size of the various 
terms. 



9. Large N and Chiral Perturbation Theory 



The large TV expansion for baryons can be combined with baryon chiral 
perturbation theory [32,59]. In the large N limit, the baryon is heavy. 
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Fig. 46. Comparison of the large A'^ predietions for the strengths of the various terms 
in the Nijmegen potential [57,58]. The ratios of couplings (relative to fp) have been 
plotted, and the expected size of the ratio is shown by a horizontal line. The shaded 
regions are an estimate of the size of the O {1/N) corrections to the leading result. 



Fig. 47. One loop correction to the pion-baryon vertex. 



and one can use the formalism of refs. [60]. The pion-baryon constant is 
order ^/N. The one-loop correction to the pion-baryon coupling constant 
is shown in fig. 47. The diagrams produce a correction to the pion-nucleon 
coupling constant of the form 

6g oc Ng^mqlogniq/ fi (9.1) 

Here gy/N is a generic pion-baryon coupling, m is a light quark mass, and 
IJ, is the scale parameter of dimensional regularization. For simplicity, I will 
denote the light quark masses by m,, and not distinguish between the u, d 
and s quark masses. One can see from eq. (9.1) that the chiral expansion 
naively breaks down in the large N limit, since the correction grows with 
A'^. This would be true if one computed the diagrams in fig. 47 by only 
including intermediate nucleon states. However, in the large N limit, the 
A-nucleon mass difference is of order and the A must be included as 
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an intermediate state as well. In the notation of section 6.1, the diagrams 
in fig. 47 arc proportional to the group-theoretic factor 

_2^»aj^jfe^ja ^ j^iaj^ifcj^ife _^ X^'' X^'' X'" . (9.2) 

The first term in cq. (9.2) is the vertex correction, and the last two terms 
are from wavefunction renormalization. The relative factor of —1/2 for the 
wavefunction diagrams arises since their contribution to the amplitude is 
\l\fZ. The sum over mesons is the sum on m, and the sum over inter- 
mediate baryons is the matrix multiplication of X. Equation (9.2) is the 
double commutator 

which is of order 1/A^^ from the consistency conditions of section 6.1. This 
converts eq. (9.1) into 

8g oc -^g^mglogmg/n, 

and the expansion parameter is now suppressed by one power of TV. In- 
cluding the entire tower of large N baryon states is crucial for this 1/N 
suppression; the double commutator is O (1) if only intermediate nucleon 
states are included. The loop expansion in the baryon sector is now a l/N 
expansion, as in the meson sector. Hadronic dynamics in the meson and 
baryon sectors becomes semiclassical in the large A'' limit [61]. The consis- 
tency conditions of section 6.1 can be derived by requiring that the chiral 
expansion has a sensible large N limit [32]. 

It is interesting to see what happens if one studies 7r-nucleon scattering 
in the large N limit without including the complete tower of large TV states. 
The nucleon is infinitely heavy in this limit, and the theory reduces to a 
strong coupling theory of pions interacting with a static nucleon, a model 
studied by Pauli and Dancoff [62]. They showed that the theory produces 
an infinite tower of states with I = J = 1/2, 3/2, . . ., precisely the spectrum 
of the baryon tower in large N QCD. The ratios of the pion- couplings are 
also the same as those of large QCD. One therefore has two possible ways 
of thinking about baryon chiral perturbation theory: (a) One starts with a 
Lagrangian with pions and nucleons. The theory is strongly coupled, with 
loop graphs of order A'^, and dynamically generates a A resonance, (b) One 
starts with a Lagrangian with pions, nucleons and A fields. The theory is 
weakly coupled, with loop graphs of order 1/A''. All predictions of (a) and 
(b) for observable quantities arc the same. The main effect of the O (N) 
strong interactions in (a) is to produce a A resonance. Once the effects of 
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this resonance are explicitly included, as in (b), the residual interactions 
are O (l/N) and weak. 

Baryon chiral perturbation theory can be formulated including interme- 
diate A states [60], and has been applied to baryon axial couplings, masses, 
magnetic moments, etc. [63-65]. Results including the A resonance pro- 
vide a good description of the experimental data. In fact, these original 
calculations showed that baryon couplings were very close to their SU{6) 
values, that there were interesting cancellations in diagrams when the A 
was included, and that ignoring the A led to disagreement with experi- 
ment. The large N approach to baryons originated in trying to understand 
these results. 

10. Conclusions 

The expansion is an extremely useful tool for understanding the prop- 
erties of mesons and baryons. S. Coleman, at the end of his 1979 Erice 
Lectures [3] remarks that "For the baryons, things are not so good. Wit- 
ten's theory is an analytical triumph but a phenomenological disaster." He 
also concludes "I feel future progress in this field rests upon constructing the 
leading approximation." Since these remarks much has been learned. The 
leading approximation (i.e. the master field) has been constructed only for 
1-1-1 dimensional chromodynamics, not for 3-1-1 dimensional chromodynam- 
ics. Nevertheless, a lot of progress has been made. Many of the applications 
to mesons phenomenology in section 4 are new. Baryon phenomenology is 
also very successful; the new results exploit large A'' spin-flavor symmetry, 
rather than the Hartree picture of Witten. 
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